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The thesis deals with the study of laser ablated carbon plasma expanding in an ambi- 
ent gas and nonuniform magnetic field using emission spectroscopy, fast photography and 
Langmuir probe. The nonuniform magnetic field and ambient gas control the plasma be- 
havior in very many ways which may be important for numerous technological applications 
and studying astrophysical processes. 

The carbon plasma was produced by focussing an Nd:YAG laser (DCR-4G, Laser 
Spectra) pulse (E = 800 mJ/pulse , A= 1.06/xm, pulse width = 8 ns) by a spherical lens on 
to the surface of a rotating graphite target kept inside the vacuum chamber. The chamber 
can be evacuated to a pressure better than 10 -4 Torr. It has a provision for introducing 
■gas into it in a controlled manner. The plasma formed was imaged on to the slit of a 
monochromator (HRS-2, Jobin Yvon) which can be tuned by a microprocessor scanning 
system. The output from the monochomator was detected by a photo multiplier tube 
(PMT)(IP28, Hamamatsu). PMT signal was recorded on the chart recorder and oscillo- 
scope (TS-8123. Iwatsu) and finally stored in the personal computer. The two dimensional 
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image of the plasma plume were recorded using intensified charged .coupled device (ICCD- 
576G/2, Princeton Instrument). The delay of the opening of the camera shutter with 
respect to ablating pulse can be varied in steps of 5 nsec. The detector consists of mi- 
crochannel plate (MCP) with a spectral response in the region 200-820 nm and 384 <576 
CCD array. The microchannel plate along with the CCD makes ICCD a highly sensitive 
device. The system essentially is an opto-electronic device where the images are stored 
electronically and analysed by the image processing software. A langmuir probe kept at a 
distance of 4 mm from the target surface at an angle of 45° monitored the evolution of the 
plasma. 

A one dimensional Lagrangian Code MEDUSA which simulates the temporal and 
spatial evolution of laser ablated plasma in vacuum is modified for the plasma expansion 
in high pressure ambient atmosphere. Plasma expanding in high pressure ambient gas 
gives rise to shock wave in the ambient gas. Rankine-Hugoniot relations for the shock 
wave are used to calculate the boundary conditions of the plasma front. The plasma 
expanding with a high velocity ~ 10 6 cm/sec acts like a hot jet streaming into an ambient 
atmosphere. The plasma boundary undergoes boundary layer mixing with ambient gas 
and gives rise to shear stress at the boundary. Due to the shear stress the plasma plume 
experiences viscous drag which decelerates the flow. Prandtl’s mixing length theory is used 
to calculate the total shear stress or the viscous drag experienced by the plasma plume. 
The external gas diffuses into the plasma plume when the pressure outside the plasma 
front is almost equal to the pressure inside the plume which gives rise to a diffusive force 
opposing the motion of the plasma. This diffusive force is proportional to the ratio of 
internal to external pressure and density of the gas. Finally, the Navier-Stokes equation 
is modified by incorporating the viscous drag term and diffusive force term. The modified 
MEDUSA (with the changed boundary conditions and modified Navier-Stokes equation) 
was used to simulate the carbon plasma expanding in high pressure (300 Torr) ambient 
gas (He, Ar, Ne and Xe). The simulated temporal variation of the plasma front for all 
the four gases matches well with the experimental observations. The simulated velocity 
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and shock wave intensity of the plasma front also match well with the experimental values. 
The plasma parameters like temperature^ pressure and electron density of the plasma front, 
can be simulated for any gas at any higher pressure with the help of modified code. 

The temporal profiles of Cl transition 2p 3 3 D° - 8f F(5/2) at 399.7 am, CII transition 
3d 2 D - 4f 2 F° at 426.7 run, CHI transition 3s 3 S - 3p 3 P° at 465.0 am and CIV transition 
at 3s 2 S - 3p 2 P° at 580.1 nm were recorded at different distances from the target surface at 
various pressures of air (0.1 mTorr. 5 mTorr and 100 mTorr). The intensity and nature of 
the temporal profile is studied for different species and for different pressure . The velocity 
of all these species was estimated from the delay of the peak of the signal with respect to the 
ablating pulse. In vacuum, it is observed that the species with higher charge have higher 
velocities ( CII ~ 1.24xl0 6 cm/sec, CIII ~ 1.76xl0 6 cm/sec, CIV ~ 1.90xl0 6 cm/sec at a 
distance of 4 mm from the target). This is because as the plasma propagates the electrons 
overtake the ions giving rise to electric field which pulls the ions resulting in higher velocity 
for higher charge. In 5 mTorr air, the velocity of CIII and CIV species decreases compared 
to that of vacuum (0.1 mTorr) value but the velocity of Cl and CII species increases. In 100 
mTorr the velocity of CIII, CIV and Cl species decreases compared to that of 5 mTorr but 
the velocity of CII species increases. The intensity of the species is observed to increase in 
presence of ambient gas. The ICCD photographs at 5 mTorr at various delay times show 
stratification of plasma plume into two components, fast and slow at 600 ns. The velocity 
of the two components is estimated to be 1.0 x 10 6 cm/sec and 7.5 x 10 5 cm/sec and 
hence termed as fast and slow components respectively. The ICCD images show that the 
fast component is bulky at earlier time (~ 800 ns) whereas slow component becomes bulky 
at later time (~ 1500 ns). In vacuum the temporal profiles show a single peak structure 
for all the species. In 5 mTorr and 100 mTorr, double peak structure is observed for 
Cl, CII and Cm transitions whereas no double peak is observed for CIV transition. The 
appearance of double peak in the temporal profiles is attributed to the stratification of the 
plasma into slow and fast components as observed in ICCD photographs. The velocity of 
fast peak increases with the increase of pressure whereas that of the slow peak decreases 
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with, pressure. At 100 mTorr, the fast peak in case of Cl and CII is stronger than that 
of the slow peak whereas in case of CIII the slow peak is stronger. At 5 mTorr. the fast 
peak is stronger than the slow peak for CL CII and CIII transitions at distances close 
to the target (~ 8 mm) whereas at later distance (~ 14 mm) it is the slow peak which 
is dominant in accordance with the results of ICCD photography. The stratification of 
the plasma plume is attributed to the diffusion of ionized air particles into the plasma 
plume and collisions with the different charged species of the plasma plume. The collided 
species give rise to slow component where as unc.oilided species give rise to fast component 
resulting in stratification. 

The ICCD photographs of the plasma expanding in a nonuniform magnetic field show 
oscillatory behavior of the plasma front at various delay times. The first contraction of the 
plasma front is seen at 150 ns and second contraction is observed at 500 ns with respect 
to the ablating pulse. The solution of energy, momentum and Ohm’s law equations of 
the plasma plume expanding in magnetic field predicts periodic bouncing of the radius of 
the plasma plume. The experimental observations are in accordance with the theoretical 
predictions though the radius of the first bounce observed experimentally (0.43 cm) is less 
than that predicted theoretically (1.64 cm). At a later time (~ 1000 ns) the plume is 
observed to break into two lobes which move toward the poles of the magnet. The ICCD 
images were also used to estimate the plasma parameters. The electron density in presence 
of magnetic field is found to be higher than that of vacuum case and shows oscillations 
because of alternate contraction and expansion of the plasmoid at 150 and 500 ns. The 
typical values of electron density at 150 ns in magnetic field, 100 mTorr air and vacuum 
are 1.0xl0 17 /cm 3 , 5.7xl0 16 /cm 3 and 2.7xl0 16 /cm 3 respectively. The temperature of the 
plasma plume is estimated assuming the plasma plume to behave as an ensemble of black 
body radiators. The typical values of temperature at 150 ns in magnetic field, 100 mTorr 
air and vacuum are 3.19 eV, 3.94 eV and 2.42 eV respectively. The temperature in case of 
magnetic field increases due to Ohmic heating caused by the surface current. In presence 
of ambient gas the temperature rises due to enhanced interparticle collisions with the gas 



particles. The calculated density and temperature of the plasma plume in case of vacuum 
are fitted with the equation of the form y = a t". The experimental observations give value 
of n = -2.85, -2.73 compared to theoretical value of n = -3, -2 for density and temperature 
respectively. 

The time integrated spectra of the plasma plume propagating in magnetic field is 
recorded at a distance of 6 mm from the target surface keeping the target at three different 
positions in magnetic field (say P, Q, R). The magnetic field strength at position P, Q, 
and R is 1.0, 3.0 and 2.7 kG respectively. When target is kept at position P, the plasma 
expands in convex magnetic field with positive field gradient, the position Q is similar 
to P but gradient is less than that in P where as in position R the field is concave with 
negative field gradient. It is found that it is not the magnetic field strength but the 
curvature and the gradient of the magnetic field which controls the plasma behavior. The 
plasma expanding in convex magnetic field with positive field gradient (position P) shows 
a maximum rise in ionic line intensity compared to other positions though the magnetic 
field is lower. Furthermore, C 2 molecular Swan band emission, d 3 7r ff - a 3 7r u , with A v = 
-1, 0, 1 is observed to dominate at this position The intensity of C 2 Swan band is very 
weak at other positions of the target. At position R where magnetic field is concave and 
gradient is negative, the intensity is less than that of no magnetic field case. Assuming 
local thermodynamic equilibrium (LTE), the temperature calculated from the ratio of line 
intensity of two lines of same ionization state at a distance of 6 mm from the target is 
highest in position P of the target (11.4 eV) and lowest in position R (4.34 eV). The increase 
in line intensity is attributed to compression of magnetic field lines ( magnetic Reynold’s 
number ~ 5) . C 2 emission in position P of the target is attributed to the movement of 
electrons due to electromagnetic interaction and that of the ions due to curvature gradient 
drift which are in the same direction favoring recombination where as they are in opposite 
direction in position R of the target. 

The time resolved spectra of carbon plasma are recorded for all the species Cl, CII, 
CIII and CIV. As the plasma expands in magnetic field, the field gradient develops at the 



plume boundary which gives rise to surface current J. The current J inside the plume is 
in a direction opposite to that at the front. The current J interacts with the magnetic 
field and gives rise to J x B force which decelerates the bulk of the plasma. But - J x 
B accelerates the boundary. Consequently, we observed three components of the plume 
viz; fast, intermediate and slow. The ICCD images show thar the plasma expanding in 
nonuniform magnetic field gives rise to two lobes, each of which has three components. 
The target is kept nearer to one pole such that two asymmetric lobes are formed and 
the three components from each lobe can be distinguished. However, for a symmetrically 
placed target between the poles it is difficult to distinguish various components. The 
observed temporal profile of the species are fitted with Gaussian functions for further 
investigations. It is found that the species at the outermost boundary, CII and CIV, show 
a single peak structure at 4 mm from the target surface whereas with increasing distance 
(~ 10 mm) multiple peaks (six peaks) are observed. They are because of two fast, two 
intermediate and two slow components arising from two asymmetric lobes. Multiple peaks 
and oscillations are observed in the temporal profile of CIII and Cl when they are at the 
outermost boundary (14 mm from the target) where CII and CIV disappear. The multiple, 
peaks and oscillations are attributed due to edge oscillations in magnetic field. 

The temporal profiles of the species show interesting features when the plasma expands 
in combined environment of ambient gas and magnetic field. The multiple peaks vanish 
and only the double peak structure is observed, similar to the one observed when plasma 
expands in only gas. Another interesting feature observed is a sudden delay in the second 
peak for CIII and CII species but not for CL This delay is not observed when the plasma 
expands in either only gas or only magnetic field. In order to investigate the delay the target 
is kept at six different positions in magnetic field and it is observed that delay occurs only 
when the plasma expands in concave magnetic field. As the target is moved further into 
the magnetic field a distinct triple peak structure is observed just before the onset distance 
of sudden delay in second peak. This behavior may be due to Rayleigh-Taylor instability 
which occurs when a plasma expands in concave magnetic field lines. The plasma (heavy 




fluid) expanding against magnetic field (light fluid) undergoes Rayleigh- Taylor instability 
when the plasma particles move along the concave magnetic fold and the centrifugal force 
drives this instability. The Rayleigh- Taylor instability along with the diffusion of air may 
be responsible for the sudden delay in second peak. The growth time of the instability is 
estimated to be 1000 ns. The effective deceleration and density gradient scale length of 
the plasma was estimated from the ICCD photographs. Since the delay is not observed 
for neutral species the possibility of recombination or charge exchange is significantly low. 

The fast ablated plasma flow into a stationary ambient medium in presence of static 
magnetic field as discussed earlier can be used to study many magnetospheric and astro- 
physical processes. At low background gas pressure (P< ITorr) collisionless phenomena 
such as those occuring when solar wind impinges on magnetospheric and ionospheric plas- 
mas and with obstacles such as comet, moons, and planets dominate. At high pressure (> 1 
Torr) collisional interstreaming plasma effects such as blast wave formation are generated 
which behave like supernova shocks. Plasma expanding in magnetic field undergo R-T 
instability. An analogous instability was observed in March 1985 AMPTE barium cloud 
experiment in the tail of magnetosphere. The plasma behavior in the curved magnetic field 
can be used to study the motion of the plasma in curved magnetic field of the earth like 
generation of ring current in the equatorial plane. Because of cosmic abundance of carbon 
in the interstellar space , spectroscopic study of carbon species is of great significance. The 
knowledge of plasma behavior in different ambient conditions elucidates our understanding 
regarding thin film deposition, surgery and many technological applications. 
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Chapter 1 
Introduction 


The importance of detailed investigation regarding laser produced plasma is increas- 
ing day by day not only beacuse of rich physics underlying the phenomena but due to its 
numerous technological applications starting from pulsed laser deposition of thin films to 
its use in providing unique laboratory test beds for investigation of many astrophysical 
processes. The latest experimental techniques like ultrafast photography and computer 
simulations using extensive computer codes have helped the scientists to reveal the com- 
plicated physics underlying the process. 

The laser produced plasma on account of its great versality that is temperature ranges 
fron 1 eV to multi keV and velocity of 10 8 cm/sec (1000 km/s) can be used as a physical 
model to study various astrophysical phenomena. The laser produced plasma flowing with 
high velocity in a static external magnetic field can be used to study astrophysical processes 
of stellar interiors. The sun emits a highly conducting plasma at supersonic speed of about 
500 km/s into the interplanetary space as a result of supersonic expansion of solar corona 
which is known as solar wind. The ablation plasma has characteristics very much like solar 
wind [1] i.e stre aming velocities of several hundred km/sec and temperatures of a few eV, 
except that the laser-plasma has much higher density and is much smaller in size than the 
solar wind. At low background gas pressure (P < 1 Torr), collisionless phenomena such as 
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those occuring when solar wind impinges on magnetospheric and ionospheric plasmas and 
with obstacles such as comet, moons and planets dominate. At high pressure (> 1 Torr) 
collisional interstreaming plasma effects such as blast wave formation [2] are generated 
which behaves like supernova shocks. Plasma (heavy fluid) expanding in magnetic field 
(light fluid) undergoes Rayleigh- Taylor instability at the boundary similar to that observed 
in AMPTE barium cloud experiment in the tail of magnetosphere (about 12 earth radii 
away) [3]. Marked striations appeared as the barium plasma expanded to its maximum 
radius. The collimation and stability properties of plasma flow across magnetic field is 
of particular relevance to the propagation of charged particle beams, solar wind evolution 
and astrophysical jets [4, 5]. These jets are similar to that observed in laboratory plasmas 
about 3 kG magnetic field via E x B drift. The propagation of laser produced plasma 
in curved magnetic field resembles the flow of charged particle in curved magnetic field 
of earth in the earth magnetosphere which gives rise to ring current flowing westward 
[6]. Understanding astrophysical phenomena such as the theory of white dwarf structure 
(mean densities ~ 10 5 — 10 8 gm/cm 3 ) requires an understanding of high density effects; 
laser plasma interaction has been utilized to measure the microscopic structure of the 
dense plasmas. Quantitative measurements of radiative opacity is of considerable interest 
to astrophysicists where laser produced plasma can be utilized [7]. Turbulence structures 
axe seen in a wide variety of space and astrophysical situations like solar wind behind 
the bow shock, supernova remnants and other astrophysical processes. These effects can 
be studied from laser produced plasma where turbulence is observed [8]. Since carbon 
has high cosmical abundance so spectroscopy of carbon plasma owes its importance to 
astrophysical processes [9]. 

Laser ablated plasma is a multifunctional fluid because of its chemical reactivity, work- 
ing gas selectivity and high energy density. Therefore, applications of plasma in material 
technology such as thin film deposition [10], production of microclusters [11], microma- 
chining of metals, plastics and ceramics [12], chemical treatment of surface layers [13], 
fabrication of microelectronic devices [14], propulsion of space vehicle [15], magnetohydro- 
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dynamic generators [16] and for environmental control have been developed extensively 
[17]. 

Laser ablation is being intensely studied and used in the computer chip industry for 
drilling printed circuit boards. Micromachining is also useful in surgery, most notably 
for kidney stones (lithotripsy) and eyes [18] (corneal and keratorefractive surgery). Laser 
ablation with very intense laser energies is being investigated as a possible X-ray source for 
X-ray lithography, as a possible commercial environmental friendly source of electricity by 
inducing fusion reactions (Inertial Confinement Fusion) [19], and as a possible laser gain 
medium for X-ray lasers. 

The application of laser ablation in the deposition of thin films onto substrates is in- 
creasingly getting popular day by day due to its various advantages over other deposition 
techniques. This technique can be used to deposit thin film of almost any material irre- 
spective of its optical properties [20] because under high power laser irradiation most of 
the materials turn optically opaque by an optical breakdown. Stoichiometry is maintmned 
between the target material and the deposited film in this process. The technique offers 
very high instantaneous growth rate and deposition can be carried out at relatively lower 
temperature due to very high kinetic energy of the laser ablated plasma species [21] . The 
process is relatively simple, inexpensive and free of contamination. It is possible to obtain 
high density films with good adhesion due to the presence of energetic species during the 
deposition process. The laser ablated deposition is mostly done in presence of gas e.g. 
oxygen for the deposition of superconducting films [22, 23] to incorporate gas in the de- 
posited films. Laser ablated carbon plasma in an ambient gas has been used for deposition 
of C6o and diamond like carbon on various substrates [24]. The dynamics of laser plasma 
is observed to change in presence of gas [25]. Thus, the knowledge of the dynamics of 
the plasma expanding in ambient atmosphere is indispensable for fabricating high quality 
films. 

Magnetic field effectively controls the local plasma thermofluid characteristics from 
both macroscopic and microscopic view points [26, 27] and hence is important from the 
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point of view of industrial applications; The "knowledge of laser produced plasma ex- 
panding in magnetic field helps to elucidate our understanding regarding beam heating 
of magnetically confined thermonuclear plasma [28], bipolar flows associated with young 
steller objects, interaction of solar winds with planetary atmosphere, formation of jets [29] 
etc. Further more, laser produced plasma in an externally applied magnetic field has been 
found to be a rich source of X-ray emission [27]. Few works have also been reported about 
change of material properties in magnetic field [30, 31]. 

It is imperative to understand the. physics of laser matter interaction in order to control 
various processes efficiently for various applications. A simplistic overview of laser matter 
interaction is presented below. 

1.1 Laser Matter Interaction 

The process of laser matter interaction is the subject of many theoretical [32-35] and ex- 
perimental [36-39] studies. The process can broadly be divided in three regimes. 


1. Interaction of laser beam with the target material resulting in evaporation of the 
surface layer (Evaporation regime) . 

2. Interaction of plasma with the incident laser beam resulting in an isothermal plasma 
and expansion (Isothermal regime) . 

3. Adiabatic expansion of the plume (Adiabatic regime). 

The first two regimes start with the laser pulse and continue until the laser pulse 
duration. The last regime starts after the laser pulse terminates. Figure 1.1 shows different 
spatial zones of the plasma along with the different temperature and density regions plotted 
on the same scale. The laser light is incident from the right hand side on to a planar solid 
target. At t < 0 the space is occupied by a cold dense material with a uniform density p s . 
For t > 0, material is heated by the laser beam and four different zones are formed. In 
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Laser 


Thermal X-rays 


Figure 1.1: A schematic diagram showing different spatial zones of the plasma during laser 
matter interaction. 

the conversion layer (CL) laser radiation is absorbed and partially converted into thermal 
x-rays. Its characteristics depend upon the wavelength of the laser light, electron heat 
conduction and the emission coefficient of the plasma for thermal radiation. It has lowest 
density and highest temperature of all the zones. It radiates equally towards the vacuum 
as well as towards the solid target. Next is the re-emission zone (RZ) where the radiation 
emitted towards the solid target is absorbed. It emits radiation through the conversion 
layer into the vacuum. Physically the re-emission occurs by multiple absorption and re- 
emission as the energy supplied by the conversion layer diffuses towards the solid target 
The re-emitted radiation is transmitted through the optically thin conversion layer and 
contributes to the thermal x-ray radiation into vacuum. It has relatively lower temperature 
and higher density in comparison to conversion layer. These two zones (I and II) are 
coupled only by radiation. Next is shock wave (SW) region where the recoil momentum 
of the expanding layer is absorbed. However, because of its large density, the energy 
tr ansmi tted to the shock wave may be neglected. Finally, the undisturbed solid which is 
solid material zone whose density is unchanged. Since there is large difference in density the 
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plasma flow can be distinguished in different zones. The details of laser target interaction 
and laser plasma interaction is described below. 

1.1.1 Laser-Target Interaction 

The interaction of laser light with the target surface depends on factors like laser fluence, 
pulse width, wavelength of the laser light and temperature dependent optical and thermo- 
physical properties of the material [20, 40-41]. For low laser irradiation (~ 10 6 Watt/cm 2 ) 
the absorbed radiation appears as heat which is evenly distributed throughout the material 
by thermal diffusion. In case of metals, the conduction electrons absorb the laser photons 
and are excited resulting in increase in kinetic energy of the electrons. Due to increase 
of kinetic energy of the electrons, the collision frequency increases which results in rise in 
temperature of the material [42]. At higher laser irradiation (> 10 6 W att / cm 2 ) , depending 
on the thermal conductivity, thermal diffusivity, reflectivity of the target material and the 
parameters of the laser pulse, intense local heating of the surface occurs resulting in a rise 
in surface temperature of the material. As the irradiation increases the temperature of the 
surface increases and a molten pool of depth ( Dt l )^ 2 is formed [20], where D is thermal 
diffusivity and 7£ is the laser pulse duration time. A further increase in the irradiation 
will cause the surface temperature of the molten pool to reach boiling point resulting in 
evaporation. This will happen when the energy deposited approximately equals the latent 
heat of sublimation [20] L, (erg/gm) ~ 1 1 t l 0 - 5 p s _1 D~ 0 ' 5 where p s is the density of the 
solid target and 1 1 is the intensity of laser irradiation. Once the vapours are formed, the 
laser light will cause further heating resulting in plasma formation. There are mainly two 
mechanisms described below which help in electron generation and plasma formation [14]. 
The first mechanism involves absorption of laser radiation by electrons when they collide 
with neutrals. If electrons gain sufficient energy they can cause ionization 

e~ + M — ► 2e~ + M + 

The electron concentration will increase exponentially with time leading to cascade break- 
down. 
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The second mechanism called a multiphoton ionization involves the simultaneous ab- 
sorption by an atom or molecule of a sufficient number of photons to cause its ionization 
or to eject an electron from the valence to the conduction band. 

M + m-phv — » M + + e~ 


If e/ is the ionization potential (or band gap) the number of photons nip must exceed 
the integer part of + 1). Multiphoton ionization is important only at short wavelength 
(A < 1 /im). 

Both cascade and multiphoton ionization require high laser irradiance usually, in ex- 
cess of 10 8 Watt /cm 2 . However, breakdown of solids has been observed at low irradiance 
~ 10 6 Watt/cm 2 in the presence of a static electric field [43]. If the solid is absorbing 
or has microscopic absorption sites, a third mechanism called thermal runaway has been 
proposed to explain the low thresholds. 

The amount of material evaporated can be calculated using energy balance consider- 
ations. The energy deposited by the laser beam on the target (El) is equal to the energy 
needed to vaporize the surface layers, the conduction losses in the target and the losses 
due to absorption of the laser energy by the expanding plasma. This is given by the energy 
threshold E t /, which represents the minimum energy above which appreciable evaporation 


occurs. The thickness of the target material evaporated/pulse (A x t ) is given by [44] 

Ax, = (l-fl)[ - - |^|| -] (cm) I 


where R, C„, AT and AH represent the reflectivity, specific heat at constant volume, 
maximum temperature rise and latent heat rise respectively. This equation is valid for 
conditions where the thermal diffusion length L t h = (2Dti) 1 ^ 2 is larger than the absorption 
length of the laser beam in the target material. 


1.1.2 Laser-Plasma Interaction 


The interaction of the laser beam with the solid target leads to vaporization of the surface 
resulting in emission of positive ions and electrons giving rise to a plasma. Once the 




1. Introduction 


8 


primary plasma is formed, the incoming laser radiation is absorbed, by a process known as 
inverse Bremsstrahlung. In this process an electron comes in the field of an ion of charge 
ze, gets accelerated, absorbs a laser photon and separates out. The absorption coefficient 
is given by 

a = 3.69 x 10 8 


z nr 


JU5 


Z/° 


1 — exp 


—hvi 


cm 


( 1 . 2 ) 


Where z, a, and T e are respectively the average charge, ion density and temperature of 
the plasma and h. kg and u L are the planch constant, Boltzmann constant and frequency 
of the laser light respectively. The term [1 - exp(-hi'£,/kgT f ,)] represents the loss due 
to stimulated emission which depends on the plasma temperature and laser frequency. 
For 1.06 fim laser wavelength, the exponential term becomes unity for T, < 13600 K 
and can be approximated by hi^/kgT., for T f . > 13600 K. The absorption term shows a 
T 0 ' 5 dependence for low temperature and T l u for high temperature. The dependence of 
frequency also changes from u L 2 to v L 2 depending on the value of (hzq,/k B T e ). Eqn.(1.2) 
shows that rate of absorption depends on nif. Tfiis results in increase in KE of the electrons 
which increases collision frequency. The energy absorbed by the electron is transfered to the 
ions by the electron ion collisions. The ions get more and more ionised resulting in increase 
in electron density. However, this process cannot go on indefinitely. The electromagnetic 
wave (laser) propagating through the plasma obeys the dispersion relation of the form 

u L 2 - cj p 2 = k 2 c 2 (1.3) 


= ( 


Att n P e 


2 1/2 


m e 


(1.4) 


Where is the angular frequency of laser light, k is the propagation vector, c, the velocity 
of light, e, the electron charge, n e , the electron density m e , the electron mass and u p is the 
plasma frequency. For ui > u p , k is real and the electromagnetic field propagates whereas 
for uji < ta p , k is imaginary and the wave is not propagated. At ui = ui p , reflection of light 
occurs at a density called the critical density n c = m e u>i 2 /Axe 2 

At the be ginni ng of the inverse Bremsstrahlung process » ui p and rate of absorp- 



1. Introduction 


9 


tion depends on n, 2 . The absorbed energy causes an increase in the kinetic energy of the 
electrons i.e. an increase in electron temperature which in turn produces further ionization 
with a consequent increase in n e and eventually approaches n c the critical density. The 
critical density extends across a plane surface some distance into the plasma. The value of 
n c is 3.15 x 10 21 /cm 3 for 1.06 fim. wavelength laser light from Nd:YAG laser. At this distance 
the plasma becomes opaque to the incoming radiation at the surface which causes laser 
light to get reflected. Laser light can no longer reach the target surface to generate new 
plasma. However, plasma growth does not cease. Because of the heating which follows the 
absorption of energy by inverse Bremsstrahlung, the plasma is driven rapidly away from 
the target surface resulting in decrease in electron density and laser light again reaches 
the surface. The energy absorbed by the electron equilibrates very rapidly so that a well 
defined electron temperature is reached. The energy gained by the elctrons is also shared 
with ions in collisions. These processes do not take place discontinuously but merge into a 
smooth self regulating regime with the generation, heating and expansion of plasma ta kin g 
place throughout the laser pulse. The expanding plasma generates pressure which drives 
a shock wave into the solid material as shown schematically in Fig.1.2. The situation is 
analogous to rocket take off condition where the hot gas exhaust result into a net upward 
thrust. 


The energy absorbed by the electron equilibrates very rapidly so that a well defined 
electron temperature is reached. The time required to attain equilibrium (the electron-ion 
thermalization time) is given by [45] 


ter 


= 2.52 X 10 8 


AT e 3/2 

n e z 2 lnA.' 


(1.5) 


Where A is the atomic mass n umb er and In A ~ 10 reprsents the coulomb logarithm for 
laser plasma. For laser produced carbon plasma, is considerably less (~ 10 -10 sec) than 
the duration of a Q-switched laser pulse ( 10 -8 sec) and so a good degree of thermalization 
is expected in the plasma. 

At higher irradiances (> 10 !2 Watt / cm 2 ) the laser ablated plasma is heated to a 
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material 


Figure 1.2: A schematic diagram showing laser plasma interaction in vacuum and rhe 
formation of shock wave moving into the target. 

very high temperature in a thin absorbing layer at the surface of the dense regime and 
has only a small optical thickness, so the self regulating regime fails to exist. Absorption 
becomes non-collisional and takes place via resonance mode conversion at critical density 7 . 
This regime has been investigated in one dimension by Fauquignan et. al. [46]. These 
plasmas are used for studying dynamics of strongly coupled plasmas [8] , short wavelength 
lasers [47, 48] and for astrophysical like plasmas [8]. 

In the initial stages of plasma expansion when the particle density is of the order of 
10 19 — 10 2 °/em 3 , the mean free path of the particles is short, and the plasma behaves as a 
continuum fluid. The equation of gas dynamics can be applied to simulate its expansion 

[49]. 

Stl 

— + '7.(nv) = 0 (1-6) 

dv 

nm{—) +- (v.V)v = —VP 
ot 


(1.7) 
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where n is the number density, m. the mass of the particle, v, the velocity and P is rhe 


pressure. The density of the plasma can be expressed as [44] 


n{x,y,z.t) = 


N T t 


v'2T' 5 nX(t)Y(t)Z(t) eXP [ 2X(t)- 2 Y{tf 2Z(tf 


( 1 . 8 ) 


where NY is the total number of ablated species at t = 17 , (pulse width). X(t), Y(t) and Z(t) 
are the dimensions of the expanding plasmas in the three orthogonal directions. Assuming 
plasma to be an ideal gas, the pressure at any point can be expressed as 

P(x*y,z<t) = n{x,y,z,t)k B T (1.9) 


It has been shown that gas dyn ami c equations follow similarity transformation [50] in 

( 1 . 10 ) 


which the velocity can be expressed as 

x dX{t), 


v = 


. , y dY(t ) , , z dZ(t ) . 

I ~T~ T ^ . \ , . J 1 , & 


X{t) dt Y(t) dt 


Z{t) dt 


The equation of motion derived from the above equations is 

X(t) 


Y dX drX 


'1 dY d 2 Y' 


'1 dZ d 2 z~ 

t dt ^ dt 2 

= Y(t) 

t dt ^ 'dt 2 

— Z(t) 

t dt ^ dt 2 




m 


t < r L . 

( 1 . 11 ) 


Where To is the isothermal temperature of the plasmas. The above equation determines 
the initial expansion of three orthogonal plasma edges. The initial dimensions of the 
plasma are of the order of mm in the transverse direction whereas in the perpendicular 
direction they are less than 1 gm. Eqn.(l.ll) shows that during initial expansion stage, 
when the velocities are small, the acceleration is very high. With an increase in expan- 
sion velocity, acceleration starts to decrease and ultimately becomes zero resulting in an 
elongated plasma shape. Since the plasma dimensions are much smaller along the target 
normal than in the transverse direction, the expansion is anisotropic in nature. 


1.1.3 Adiabatic Expansion of the Plume 


After the termination of laser pulse, no particles ared evaporated or injected into the inner 
edge of the plasma and the expansion is mostly adiabatic. The thermal energy rapidly 
gets converted into kinetic energy giving higher expansion velocity of the plasma. Using 
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the adiabatic equation of state for the expansion of plume after the termination of laser 
pulse, we have 

r[X(t)y(t)Z(t)] 7_1 = Constant (1.12) 

1 _ 

P 

_ .u.VT = (1 - 7 )rv.u (1.14) 


dF_ 

dt 


+ V.VP 


2 

n 


dn 

— + v.Vn 


= 0 


(1.13) 


where T is the temperature, n and P remains similar as before 


n(x, y, z, t) = 


N t 


V27r L5 X(t)Y(t)Z(t) 


exp 


r 


2 x(ty 2 Y{ty 2Z{ty\ 


(1.15) 


Solving the above equation one obtains the following equation of motion for the adiabatic 
regime, 


X(t) 


'drX~ 


'<PY' 


~d z Z' 

kT n ' 

. dt ' 2 . 

= y(t) 

dt 2 

= Z(t) 

dt 2 

m. 


X 0 Y 0 Z o 


17-1 


[X(t)Y(t)Z(t)\ 


t > t l . (1.16) 


where Xq, Y 0 , Z 0 are the initial orthogonal edges of the plasma after the termination of 
laser pulse (t =77 ) that is after the isothermal regime. In this regime, the initial transverse 
dimensions are much larger than rhe perpendicular dimension which correspond to the 
expansion length in isothermal region. As the velocities are dictated by these lengths, the 
highest velocity is in the direction of smallest dimension. This gives rise to the character- 
istic plasma shape elongated outward from the surface. 


1.2 Computer Simulation 

The description of laser matter interaction given in the last section is very simplistic 
and does not take into account many important processes like loss of energy through 
Brehsstrahlung radiation by electrons etc. The studies related to laser produced plasma 
involve variety of complex phenomena like driver (laser or ion beam) interaction with 
the plasma, charged particle transport and collisional processes, nuclear fusion reaction 
kinematics, hydrodynamics and thermodynamics of the plasma etc. The simultaneous 
presence of these phenomena neccessitates invoking of numerical methods in the form of 
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computer codes which simulate various processes taking place in the plasma. Each and 
every process having some significance is considered in the model and is dealt with one by 
one along with their inter-effects. The term simulation is referred to as "It is an activity 
whereby one can draw conclusions about the behavior of a given system by studying 
the behavior of a corresponding model whose cause and effect relationship is same as 
(or similar) to those of the original system". Computer simulation study is actually a 
form of numerical experimentation in which mathematical model takes the place of an 
experimental facility. In the past, several computer codes have been developed with the 
ultimate aim of studying laser produced plasma and inertial confinement fusion. The name 
of few of them are MEDUSA [51], CASTOR [52], LASNEX [53], MINIRA [54], ILESTA 
[55], HEMICO [56], MULTI (1-D) [57] etc. We have modified 1-D Lagrangian code Medusa 
(version MED-101) for our work [58]. 

1.3 Laser Ablated Plasmas in the Presence of an Am- 
bient Atmosphere 

The interaction of laser ablated plumes with background gases is receiving increased at- 
tention due to its increasing need in technological applications and is the subject matter 
of theoretical and experimental investigations. Knight [33] developed a theoretical model 
where homentropic expansion of vapor against ah ambient atmosphere has been examined 
applying conservation of mass, momentum and energy across the Knudsen layer. Bel- 
lantone et a 1 [59] generalised the model where entropy can vary with time that is the 
homentropic assumption of previous case is relaxed and the more general case of flow can 
be examined but its usage is restricted to a relatively low intensity regime. Aden et al 
[32] extended the previous case from one dimensional to three dimensional model which 
matched well with their experimental data and showed that there exists a material depen- 
dent minim al laser intensity where the vapor properties at the vaporization front become 
independent from further expansion of the plasma against the ambient gas. Macfarlane 
et al [34] had given a one dimensional model where radiation hydrodynamics calculations 
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are performed to study the structure and evolution of blast waves resulting from laser - 
generated plasma expansions for the ambient pressure ranging from 25 mTorr to 5 Torr. 
Giuliani et a 1 [35] had given a numerical model where laser target and debris-background 
interactions are modeled solving mass continuity, total momentum and separate ion and 
electron internal energy equations. The model is appropriate for background pressures 
greater than 1 Torr. Ananin et ai [60] had given a model for plasma expansion at low pres- 
sure ambient gas of the order of 0 - 1 Torr. The analysis is based on multifluid gas dynamics 
approach because of diffusion of external gas particles into the laser plasma. The blast 
wave formation and shock is well described analytically by self similar solutions derived 
by Sedov [61] which works well at higher pressure. The model is true for an instantaneous 
point explosion and with negligible radiation loses. However, the model is insufficient for 
getting the detailed thermodynamic and hydrodynamic properties of the blast wave. In 
low gas pressure the plume propagation could be described by Monte Carlo simulation 
[62]. Recently a numerical model has been proposed where Navier Stokes equation for the 
case of vacuum is modified by incorporating viscous drag force and diffusive force. The 
boundary conditions relate to Rankine-Hugoniot conditions for shock wave [58]. 

The behavior of laser produced carbon plasma expanding in an ambient atmosphere 
has been investigated experimentally in the recent past. Carbon plasma expanding in an 
ambient air [39], helium [63], argon [64] atmosphere has been studied extensively. Line 
emission enhancement in laser-ablated carbon plasmas in ambient helium and argon gas 
at -various pressure is studied by Abhilasha and Thareja [65]. Role of ambient gas (helium 
and argon) on laser ablated plumes for thin carbon film deposition has been discussed 
by Thareja et a 1 [66]. Geohegan et ai [67] has discussed the role of He, Ne, Ar and Xe 
at high pressure (300 Torr) on the propagation dynamics of laser ablated graphite. The 
interaction of laser ablated plumes with background gas is receiving increased attention 
due to its importance for film growth by pulsed laser deposition [68-75] particularly for 
deposition of high T c superconductors [76-78], polymeric thin films [79], diamond like 
carbon films [80-82], TiN films in presence of N 2 [83] and ceramic thin films [84]. At higher 
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gas pressures laser ablation can be utilized to produce novel clusters such as fullerenes 
and doped fullerenes [85, 86]. Real time investigation of the plume dynamics during 
laser ablation of germanium in an oxygen environment has been reported by Vega et ai 
[87]. Correlation between target-substrate distance and oxygen pressure in pulsed laser 
deposition of YBa 2 Cu 3 07 has been discussed by Kim et al [88]. Anomalous enrichment of 
C 2 ~ ions by the laser ablation of graphite in Ar jet has been observed by Park et al [89]. 
Formation of shock wave structure during expansion of a laser plasma in a low-density 
[90-92] and high density gas [93] has also been investigated. 

There have been several reports on propagation of laser produced plasmas in ambient 
atmosphere, however, the dynamics of carbon plasma expansion in ambient gas at low 
pressure is still not clear. The principal diagnostics technique like emission spectroscopy 
[34, 94-97], laser induced fluorescence [38, 98], langmuir probe [99-102], streak photography 
[103, 71], ICCD photography [104, 105] etc have been used to study laser produced plasmas. 
The plasma propagating in low pressure. ambierit atmosphere splits into two components, 
fast and slow components [25]. Some researchers proclaim that Rayleigh-Taylor instabil- 
ity occuring due to the difference in density of the propagating plasma and the external 
ambient gas to be the major cause of this stratification into slow and fast components 
[103, 25]. For high background pressure and at later time a blast wave model describes 
well the position of luminous shock front whereas at low pressure and at early time better 
fit is obtained from classical drag force model [78]. However, these models fail to predict 
the fast component of the propagating plume. Attempts have been made to simulate the 
occurence of slow and fast components numerically by introducing a dynamic source effect 
that accelerates the expansion of laser ablated material in the direction perpendicular to 
the target [106]. A theoretical model based on a combination of multiple inelastic scatter- 
ing and hydrodynamic formulation gives good fit to experimental data for Si in He and 
Ar [100]. The stratification of the plasma plume as observed by fast photography is in 
accordance with time resolved emission spectroscopy [105]. 

The plasmas produced in an ambient atmosphere differ significantly from that pro- 
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duced in vacuum. Laser plasma expansion in low density gas can broadly be divided into 
two categories viz. 0-1 Torr and above 1 Torr. At low pressure (0-1 Torr), the mean free 
path of plasma ions either exceeds or is comparable with the dimension of the experimental 
chamb er (~ 10 cm), there is a strong interpenetration of the laser plasma and the ambient 
low density gas [60]. Above 1 Torr pressure the laser plasma expansion results in the ap- 
pearance of a shock wave structure in the gas [106]. Various processes like recombination, 
charge transfer, collisional and collisionless interactions, shock/blast wave interactions, mi- 
croinstabilities in the laser plasma occur during laser plasma expanding in ambient gas. 
They are briefly described below. 

1.3.1 Recombination Process 

When an ion and electron collides, particularly at low relative velocity, they have a finite 
probability of recombining into a neutral atom. To conserve momentum a third body 
must be present. If this third body is an emitted photon, the process is called radiative 
recombination. If it is a particle, the process is called three-body recombination [107, 108]. 
For the radiative recombination 

C 7 ' + e = C (9 ~ 1)+ + hv 


The radiative recombination rate being 


given by 
= -a R n e m 

R 


where a R is the radiative recombination coefficient, 

a R = 2.7 x 10 ~ l3 z 2 T e ~ 3/4 cm 3 /s 


(1.17) 


( 1 . 18 ) 


where 14 and ip are the electron and ion densities, respectively: T e is electron temperature 
in eV and z is ionic charge. Three body recombination process is defined as 


C q ~ + e + e = C (? ~ 1)+ + e 
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(1-19) 

cm 6 /s (1.20) 

From Eqs. (1.18) and (1-20) it follows that radiative process is important only close to the 
target surface and the three body recombination is a dominant process beyond a few mm 
from the target surface. 

1.3.2 Charge Transfer 

When an ion enters into the vicinity of a neutral/ionic particle, there is a finite possibility 
that an electron may be transferred from the neutral/ionic particle to the other ion and 
thereby exchange takes place [109]. This process is referred to as charge transfer and is 
represented by 

X~ +Y X + Y + ±AE (1.21) 


Three body recombination rate R, being given by 


Rc 


dn e 


dt 


— —accUr'-n^ 


where a c is the three body recombination coefficient, 

a c = 9.2 x lO-' 27 z 3 TY V2 ln[(z 2 + 1) 1/2 ] 


The quantity A E is called the energy defect in the electronic collision. The energy de- 
pendence of the cross section for a particular charge-transfer reaction depends upon the 
magnitude of the energy defect. The charge-exchange rate of ions is given by the relation 

R = rii < Vicr > (1-22) 


where V; is ion streaming velocity, and cr the cross-section for the charge transfer. The 

(1.23) 


cross section for the ion-ion charge transfer interaction is given by [109] 

A-KZi 2 zo 2 eHnA. 

a{zx,z 2 ) = 


nWil + ma/mii 


where index 1 refers to the laser plasma ions, 2 to the ambient gas ions, z is the charge 
number, In A ~ 10 is the Coulomb logarithm for the laser plasma, v the mean velocity of 
laser plasma ions, nii the mass of ion and /i the reduced mass of ions. The cross-section cr, 
for the conditions of our experiment is ~ 10 -18 cm 2 . Seely et. a l. [110] have examined the 
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charge exchange process occurring between C VTI and Ar III in expanding laser produced 
carbon plasma in details. 

1.3.3 Shock / Blast wave 

In case of plasma expanding in an ambient gas (~ 1 Torr), the laser plasma acts like a 
pl ung er which forms a shock wave in the gas and there is a contact boundary between rhe 
gas and the laser plasma as shown in Fig. 1.3. The background gas is compressed into a 
thin shell at the shock front. As the shock front expands, more and more of the background 
gas is swept up by the shock wave resulting in decrease in the expansion velocity of rhe 
laser produced plasma. This kind of unsteady shock wave is known as blast wave. For 
high background pressure a blast wave model [50, 111, 112] has been used to describe rhe 
luminous shock front caused by the expansion of laser produced plasmas. This model is 
s im ilar to a shock wave caused by an point explosion with an explosive release of energy Ei, 
through a background gas (density po). The propagation of shock front by the background 
gas follows the distance-time relationship 

St = C 0 (E L /p 0 ) 1/5 ^ 5 (1.24) 

where Co is a constant. The above expression for the shock front holds good for the 
expansion in high pressure ambient gas. At low pressure a classical drag model shows 
better agreement where the ejected pulse of ablation products regarded as an ensemble 
experience a viscous force proportional to its velocity through the background gas. The 
equation of motion is 


1 

II 

<3 

(1.25) 

v = v 0 exp(-/3 s t) 

(1.26) 

x - x f [l - exp(-/3,t )] 

(1.27) 

where 3 S is the slowing coefficient. v 0 is the initial velocity and xy 

distance of the plume. 

= v 0 / (3 s is the stopping 
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Figure 1.3: A schematic diagram showing the expansion, of laser produced plasma and rhe 
formation of shock wave into the ambient gas. 

1.3.4 Microinstabilities 

Due to difference in density of the two immisible fluids the interface develops some kind of 
corrugations or perturbations which if grow with time give rise to instabilities. There are 
three standard methods normally used to determine the stability or the properties of the 
instabilities associated with a given plasma equilibrium [113]. 

The simplest is the intuitive approach, in which the equilibrium is subjected to a per- 
turbation that alters the forces acting on the plasma. If these modified forces act to increase 
the initial perturbation the plasma is unstable. This approach does not give quantitative 
information about growth rates, but it does give insight into the instability mechanism and 
show which equilibrium situations require further analysis, as well as providing information 
about which modes are likely to be unstable. 

The second approach is based on energy considerations. An equilibrium state of a 
plasma is stable if potential energy is a minimum and unstable otherwise. To investi- 
gate the stability of a plasma with this method it is necessary to calculate the change 
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in potential energy of the plasma due to a given perturbation. The plasma equilibrium 
is stable if the change in potential enegy is positive for all perturbations allowed by the 
plasma equations. In the presence of even one allowed perturbation for which the change 
in potential energy becomes negative, the plasma becomes unstable. Using this criteria it 
is possible to estimate the values of parameters of the plasma which take the system ro 
unstable state. However, it is not a useful method for determining the growth rate of a 
particular instability. 

A third approach to determine the stability of a plasma equilibrium is normal mode 
analysis. It assumes that if the equilibrium under investigation is perturbed, the linearized 
plasma equations for the time development of the perturbation can be solved, subject ro 
the appropriate boundary conditions, assuming a time dependence of the form exp(-i -j, 
t) . This procedure gives an equation for u\ in terms of the equilibrium parameters. The u;, 
’s from this equation may be real, imaginary, or complex. If ail the tO\ ’s are real, then the 
perturbation variables oscillate harmonically and the plasma is said to be in a stable state. 
If any or all of the uq ’s have positive imaginary parts, the system is unstable. However, 
the normal inodes will grow in time. Normal-mode analysis provides complete information 
about the instabilities associated with a particular plasma equilibrium. The development 
of any initial perturbation can be followed up to the limits imposed by the linearization 
of the equations. Unfortunately, normal-mode analysis can be applied only in those cases 
where the plasma equilibrium is simple enough to allow solution of the plasma differential 
equations. 

One common example which often occurs at the interface of two moving fluids is 
Rayleigh Taylor (RT) instability. RT instability occurs at the interface of the two fluids 
of different densities when accelerated in a perpendicular direction to their interface; this 
interface is unstable or stable according to whether the acceleration is directed from heavier 
to the lighter medium or vice-versa [114, 25]. In other words Rayleigh Taylor instability 
occurs when the density grows in the direction against the accelerating force. A heavy 
fluid supported against gravity by a light fluid is a common example of Rayleigh Taylor 
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instability. The details of the Rayleigh Taylor instability and the calculation of growth 
using the third approach is presented in Chapter 6. 

1.4 Equilibrium in Plasmas 

In complete thermodynamic equilibrium [45] (with both kinetic and radiative equilibrium 
at the same temperature T) every atomic processes occuring in the plasma is balanced by 
an equal and opposite process. In this case the following conditions are satisfied : (a) all 
particles, electrons, neutral species and ions obey the Maxwell velocity distribution law; (b) 
the population distributions over the states of any atom or ion are given by the Boltzmann 
formula; (c) the number of ions in stage z relative to the number in stage (z-1) is given by 
the Saha equation, (d) the intensity distribution of the radiation as a function of frequency 
and temperature is given by the Planck’s formula. All these conditions are characterized by 
the same temperature T. Rarely, especially in terrestrial plasmas, a perfect thermodynamic 
equilibrium is approached; the very fact that radiation is emitted from a plasma prevents 
thermodynamic equilibrium from occuring. 

The commonest plasma model conforms to local thermo dynamic equilibrium (LTE). 
In LTE it is assumed that collisional events, in particular the events involving electrons, 
determine the behavior of the system. The collisions are assumed to be governed by the 
same laws that hold in total thermodynamic equilibrium. The temperature used is the 
electron temperature, T e , since electrons dominate the collision process. The radiation 
distribution however is not given by Planck’s equation and radiative effects are assumed to 
play an insignificant role in determining the equilibrium in the plasma. For LTE to hold, 
the elctron density must be sufficiently high. According to McWhirter [115], for LTE to 
hold lie, the number of electrons per cm 3 must satisfy; 

n e > 1.6 x 10 12 T e x/ V (1-28) 

where x is the excitation potential (in eV) of the transition. At low electron densities a 
second limi ting type of equilibrium called coronal equilibrium (so called because it describes 
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conditions in the solar corona) may be attained. In this case excitation and ionization occur 
as the result of electron collisions whereas de-excitation and recombination occur bv the 
e mis sion of radiation. The plasma is assumed to be optically thin so that re-absorprion of 
the emitted radiation is negligible. 

I. 5 Laser Plasma in Magnetic Field 

The physics of plasma flow across magnetic field lines is important in many laboratory and 
space plasmas. The familiar qualitative picture of a plasma moving in transverse magnetic 
field fall into one of the two extreme categories 

(1) In case of a low pressure plasma, the magnetic field diffuses into the plasma resulting in 
polarization of the plasma, i. e. an electric field E is established which is perpendicular to 
both B and the plasma velocity. The plasma moves across the B field by E x B drift. This 
case has been discussed qualitatively by Schmidt [116] and a more complete discussion has 
been given by Bostick [117]. 

(2) In case of high pressure plasma, the B field is excluded from the plasma interior. If 
the plasma pressure exceeds the magnetic pressure, the plasma "punches" its way across. 
To produce a field-free cavity of volume V in a field B requires energy W which is scored 
as potential energy in the plume, given by 

W = f-V (1.29) 

2/i 0 

where n o is the permeability of free space and B is expressed in kG, V in cm 3 and W in 

J, W = 4 xlO -3 B 2 V. The maximum volume to which the plasma can expand is found by 
setting W equal to the plasma energy [118]. As the plasma expands in the field, currents 
are generated in the surface which result in O hmi c heating. In this way some of the 
expansion energy is converted into thermal energy. The interface between the plasma and 
the region of greater field strength is unstable against small perturbations, another form 
of Rayleigh-Taylor instability. 

The expansion of laser produced plasma expanding in magnetic field has been dis- 
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cussed theoretically by several authors. The simple model of a cylindrical plasma which 
expands freely, parallel to its axis, along the magnetic field lines was discussed by Shkuropat 
and Shneerson [119] whereas that across a magnetic field was investigated by Bhadra [120] 
who for simplicity assumed a spherically symmetric uniform magnetic field. Bernstein and 
Fader [121] considered the case of a spherical, resistive, collision dominated plasma, which 
was assumed to retain its spherical shape as it expands in a uniform magnetic field. Haught 
et. al. [122] included a viscous term and calculated the temperature, position and velocity 
of the plasma boundary from a particle of LiH, 40 /xm in diameter, in a magnetic field of 
10 kG. The expansion of a uniform spherical plasma into a uniform background plasma in 
a uniform magnetic field was discussed by Wright [123] who found an analytic solution for 
the magnetic field configuration that is how the magnetic field lines change during the early 
stages of the expansion. Hains [124] showed that departures from azimuthal s ymm etry ran 
impact a rotation to an isolated plasma in a magnetic field. 

Cavaliere et. al. [125] suggested that a conducting plasma created inside a fixed 
conducting cylinder compresses the field lines if it expands perpendicular to the field lines. 
This configuration increases the effective field and has a strong confining effect on the 
plasma. Martineau and Tonon [126] investigated the developments of a plasma heated by 
a laser pulse of 6.6 GWatt peak power in a magnetic field of 1 MG within a conducting 
cylinder of radius 2 cm. Oscillations were observed along the minor axis, perpendicular to 
the field whereas plasma expanded along the major axis, along the field. The contribution 
of inverse Bremsstrahlung and of induced Compton effect to the absorption of energy from 
the laser beam were both calculated and found to be roughly equal. 

Schwirzke and Tuckfield [127] have pointed out that if a plasma expands sufficienrly 
in a magnetic field for the collision frequency to be small, two-stream instabilities may 
develop which increase the resistivity of the plasma and thus reduce the rate of cooling 
by thermalizing some of the energy of expansion. Dunn and Lubin [128] discussed the 
radiation emitted from a plasma expanding in a magnetic field. In addition to normal 
cyclotron radiation from the volume of the plasma, whose intensity varies as B 2 and as 
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T e , a surface contribution was calculated whose intensity varies ;is B and as Lubin 

suggested that the time history and the extent of the interaction of the expanding plasma 
could be determined from a study of the emission spectrum in the vicinity of cyclotron 
frequency. Behavior of plasma moving in nonuniform magnetic field showed that the 
plasma plume splits into two lobes which move towards the poles of the magnet [1291. 

In eighties and nineties extensive work had been reported on plasma instabilities in 
magnetic field and plasma jets. The plasma expanding in high magnetic field (> 11 kG) 
gives rise to plasma jets. Mostivych and Stamper [130] had observed highly collimated 
plasma jets by laser irradiation of solid barium targets. The plasma jet was observed to 
narrow or "focus" in the plane perpendicular to the field, while in the plane of the field 
the plasma expanded along the field lines and displayed fiutelike striations. The narrowing 
of the plasma jet was explained in terms of the configuration of the plasma polarization 
fields, while flute structure was identified as an electron-ion hybrid (E.I.H) velocity-shear 
instability. E.I.H instability was also observed by Peyser et. ai. [131] where intense 
laser pulse was focussed onto the interior wall of a hollow glass cylinder to provide some 
degree of collimation to the emerging plasma. Fast framing camera photographs showed 
markedly different geometrical features for plasma expanding in magnetic fields from B 
= 0 kG to B = 10 kG. At 10 kG strong jets were observed whereas below 3 kG no jet 
was observed. Dimonte and Wiley [132] studied the expansion of laser ablated plasma 
in a magnetic field with a Faraday-rotation magnetic imaging probe and fourier-analyzed 
optical plasma images and observed plasma instabilities which evolve from short to long 
wavelengths. Nishiyama et. ai. [133] have studied the behavior of an argon plasma jet 
in external electric and magnetic field experimentally and through numerical simulation 
using three fluid model. Ripin et. ai. [134] has observed linear and nonlinear features 
of a strong plasma-magnetic field interchange Rayleigh-Taylor instability in the limit of 
large ion Larmor radius, which undergoes rapid linear growth culminating in face-shearing 
flute tips. Okada et. ai. [135] observed flute-like instability in a cylindrical column of 
laser produced plasma generated in a uniform magnetic field. Kent et. a 1. [136] observed 
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transverse Kelvin- Helmholtz instability which occurs in a fluid with a gradient in mass 
flow velocity, in a rotating plasmas (Q machins). Ganguli et. al. [137-141] examined 
in detail the instability of a plasma characterized by a sheared flow transverse to the 
ambient magnetic field. Lehnert [142! discussed the boundary conditions of a magnetically 
confined plasma and related induced surface current which are generated by free-boundary 
displacement currents across an inhomogeneous externally applied magnetic field. The 
study is important from the point of stability criteria. Recently signature of Rayleigh- 
Taylor instability of a plasma expanding in concave magnetic field has been observed 
[143]. 

A few reports have been made regarding material studies in magnetic field. Dirn- 
berger et. al. [30] observed enhanced excitation and ionization of magnesium plasma in 
presence of magnetic field. The results are interpreted in terms of a magnetohydrodynamic 
interaction. Lash et a I [31] had observed enhancement in the ionization due to increased 
electron confinement [144] in a laser ablated aluminium plasma in a transverse magnetic 
field. Significant changes in optical emission spectra were observed in laser ablated carbon 
plasma expanding in magnetic field [70]; enhanced collisions between ablated species due 
to cyclotron motion are held responsible for the observed change in behavior. Surprisingly, 
all these experiments were done in nonuniform magnetic field but the effect of nonunifor- 
mity of magnetic field was not investigated. Recently we have observed that curvature and 
gradient of magnetic field greatly controls the plasma behavior [145]. 

1.6 Diagnostics 

In the present work, laser produced carbon plasma expanding in different ambient condi- 
tions had been studied by emission spectroscopy, ICCD photography and Langmuir probe. 
Since many different diagnostics have been used to study laser ablated plasma, a summary 
of the diagnostics utilized and the parameter that each measures are presented in Table 
1 . 1 . 
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Table 1.1: Various diagnostics and the corresponding plasma parameters. 


Diagnostics 

Plasma parameters 

Emission Spectroscopy 
[94-97, 146-153] 

Electron temperature 

Electron density 

Ion Probe Diagnostics 
[101, 102, 154-159] 

Electron temperature 

Electron density 

Photography k Imaging 
[71, 103-105, 175] 

Velocity of the plasma plume. 
Electron temperature 

Electron density 

Laser Induced Fluorescence 
[99-100, 170-174] 

Electron temperature 

Time of Flight Mass Spectroscopy 
[160-163] 

Velocity of the species, 
states of ionizations 

Absorption Spectroscopy 
[164-169] 

Electron density , temperature 

Interferometry 
[130, 177] 

Electron density 

Laser beam Deflection Method 
[179] 

Density gradient 
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1.7 The Present Work 

In the present work the dynamics of laser ablated carbon plasma expanding in vacuum, 
low pressure ambient atmosphere and nonuniform magnetic field has been reported. The 
diagnostics used to reveal the behavior of laser produced carbon plasma in different am- 
bient conditions are emission spectroscopy - the time integrated and time resolved, ICCD 
photography and Langmuir probe. 

Chapter 2 deals with the details of the experimental techniques and instruments used 
in the experiments. 

Chapter 3 deals with computer simulation studies of laser produced carbon plasma 
expanding in high pressure ambient gases. One dimensional Lagrangian Code Medusa 
which simulates the temporal and spatial evolution of laser ablated plasma in vacuum is 
modified for the plasma expansion in high pressure ambient atmosphere [58]. The temporal 
evolution of simulated plasma front for the case of plasma expanding in He, Ne, Ar and Xe 
gas at 300 Torr matches well with the experiment. The modified code can also simulate 
the plasma parameters viz. pressure, temperature, density of the plasma expanding in any 
gas at any higher pressure. 

The bahavior of laser produced carbon plasma expanding in vacuum and low pressure 
ambient atmosphere [105] (5 mTorr and 100 mTorr) is presented in Chapter 4. The tem- 
poral profile of Cl transition 2p 3 3 D° - 8f F(5/2) at 399.7 nm, CII transition 3d 2 D - 4f 2 F° 
at 426.7 nm, CUE transition 3s 3 S - 3p 3 P° at 465.0 nm and CIV transition 3s 2 S - 3p 2 P° at 
580.1 nm were recorded. The transitions show a single peak structure in vacuum. In low 
pressure ambient atmosphere they show a double peak structure except that in CIV. The 
nature of double peak structure is different for different species and changes from 5 mTorr 
to 100 mTorr ambient pressure. The ICCD images of plasma plume show stratification of 
the plasma plume at 5 mTorr. Based on the above observations the formation of double 
peak structure and stratification of the plasma plume is explained. 

Chapter 5 describes the behavior of laser produced carbon plasma in nonuniform 
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magnetic field in vacuum. The ICCD photographs of the plasma expanding in a nonuniform 
magnetic field show oscillatory behavior of the plasma front at early time as predicted 
theoretically [ITS]. At a later time (~ 1000 ns) the plume is observed to split into two lobes 
which move toward the poles of the magnet. The ICCD images of the plasma plume are 
used to calculate the plasma parameters viz. electron density and temperature as a function 
of time for plasma expanding in vacuum, ambient gas and magnetic field. Oscillations are 
observed in plasma parameters. Multiple peaks are observed in the temporal profile of Cl. 
CII, CIII and CIV transitions because of formation of two asymmetric lobes [129]. Each 
lobe has three components viz. fast, intermediate and slow. The time integrated spectra 
of carbon plasma is taken keeping the target at different positions in magnetic field. It is 
observed that the curvature and gradient of magnetic field greatly affect the behavior of 
emission spectra [145]. The emission lines of ionic species as well as neutrals are excited 
and molecular C > Swan band emission is observed when the plasma moves in convex part of 
magnetic field with positive field gradient. The observed change in emission characteristic 
is attributed to the compression of the field lines due to moving plasma and recombination 
of electrons due to curvature gradient drift. 

The expansion of laser produced carbon plasma in nonuniform magnetic field in pres- 
ence of ambient atmosphere is presented in Chapter 6. A delayed second peak is observed 
which is not present in case of plasma expansion in only ear or only magnetic field [143]. By 
keeping the target positions at different regions of magnetic field it is observed that delayed 
second peak occurs only in that region of magnetic field where field lines are concave to 
the moving plasma. The possibility of Rayleigh Taylor instability along with diffusion of 
external air as the cause of delayed second peak is discussed. 

Chapter 7 s umm arizes the results of present work. 




Chapter 2 

Experimental Techniques 


An overview of various experimental techniques used to study laser ablated plasma 
is described in Chapter 1. Various experimental techniques used to study the dynamics 
of laser ablated plasma expanding in vacuum, low pressure ambient gas and nonuniform 
magnetic field are described in detail in the present section. 

Nd.YAG (Spectra Physics, Model DCR-4G) laser along with its harmonics was used 
for laser ablation. It delivers energy upto 1 mJ/pulse with full width at half maxima 
(FWHM) of 8 ns at fundamental with a repetition rate 10 pulses per second. The laser 
beam has a Gaussian limited mode structure, with beam divergence less than 0.5 mrad and 
spot diameter of 8 mm. The energy of the laser was monitored using a laser power meter 
(Ophir Model 30A) by placing the meter in the path of the main beam. Laser energy was 
varied by changing the voltage on the laser oscillator and amplifier. Burn patterns taken 
at different energies show no significant variation in the mode pattern. 

The interaction chamber and the experimental configuration used for ablation studies 
is shown in Fig.2.1. The laser beam was focused onto the spectroscopic grade pure graphite 
target rod attached to the target holder in the chamber. When laser ablation studies were 
done in presence of magnetic field, the target was kept in the magnetic field as shown 
in the inset of Fig.2.1. The target rod was continuously rotated and translated with the 
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PRINTER 


Figure 2.1: Schematic of the experimental set up used in the experiment. Position of 
the magnet and target is shown separately. T: Target, MC: Monochromator, SO: Storage 
Oscilloscope SCR: Strip Chart Recorder, PG200: Pulse Generator, L1/L2: Convex Lens, 
ICCD: Intensified Charged Coupled Device. LP : Langmuir Probe. The magnetic field 
along Z and X axis is shown in the inset. 
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help of an external stepper motor so as to ailow each laser pulse to. interact with the fresh 
target surface. The interaction chamber is made up of mild steel with eight output ports. 
One of the output ports is used for laser entrance and the two perpendicular ports were 
used for recording the optical emission from the ablated plume. The interaction chamber 
was evacuated to a pressure better than ID -3 Torr using a rotary / diffusion pump. The 
ambient gas was fed into the chamber through a needle valve connected to the interaction 
chamber. Experiments were carried out in presence of air at various pressures ranging 
from 10~ 3 - 10 -1 Torr. 

2.1 Optical Emission Diagnostics 

A high temperature and high density plasma was created by focussing the laser beam 
using a 50 cm lens on to the pure graphite target. The plasma radiation emitted from 
the target was imaged onto the entrance slit of a monochromator (Jobin Yvon, HRS- 
2) with a lens, so as to have one to one correspondence with the plasma and its image 
onto the slit of the monochromator. The monochromator was continuously tuned using 
a microprocessor controlled scan system [180] . The output from the monochromator was 
detected with a photomultiplier tube (Hamamatsu, IP28) and recorded on a strip chart 
recorder or displayed on the screen of the storage oscilloscope (Iwatsu, TS-8123) interfaced 
with PC. The recorded optical emission data are used to evaluate the electron temperature 
using the emission lines intensity ratio, electron density using Stark broadened profile of 
a transition and the vibrational temperature using the intensity of Swan band heads. 
Investigation of the optical emission of the palsma plume also gives information about the 
optical and temporal evolution of the species produced such as excited atoms, ions and 
molecules during laser target interaction. 

2.1.1 Electron Temperature 

Various methods of calculating electron temperature using spectroscopic observations are 
described in the literature [181, 182]. We have used the relative intensities of the spectral 
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lines from a given ionization state for estimating electron temperature. Assuming plasma 
in local thermodynamic equilibrium (LTE), the electron temperature can be estimated 
using, [181, 183] 


keT e = 


E' - E" 

hi! W) 


( 2 . 1 ) 


where E" and E' are the excitation energies of the lower and upper levels, respectively. 
I" and I' s are the intensities, g" and g' are the statistical weights, and A" and A' are the 
transition probabilities and A" and A' are the wavelengths of various transitions, k# and 
T P represent the Boltzmann constant and electron temperature, respectively. The slope of 


the curve (E' - E") against j^Zn 


gives electron temperature, the parameters used 
for estimating electron temperature are available in the literature [184] . The measured 
electron temperature using intensity of CIII species in our work lies between 4-7 eV. 

2.1.2 Electron Density 

The Stark broadened profile of a transition for which the Stark-broadening parameters 
are available in the literature, can provide electron densities of the plasma [39. 185, 186]. 
Stark-broadening arises due to the coulomb interaction of the emitted species with both 
the electrons and ions expanding in the plume. The line profiles were recorded keeping the 
monochromator resolution to maximum [40]. Observed lines with true half width being 
given by, [182, 36] 

AA !r .^ AA observed AXinstrument (2-2) 

where AA of a line (FWHM) is given by 


AA = 2W ‘ + 3 ' 5 ^(^) ' X (1 - (^i?) & ( 2 - 3 ) 


First term in the Eqn.(2.3) gives contribution from electron broadening and second term 
is the ion broadening correction; W e is the electron impact parameter which can be inter- 
polated at different temperatures and A, is the ion broadening parameter, both W e and Ai 
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are weak functions of temperature, m is the electron density and Np is number of particles 
in Debye sphere given by 

N d = 1.72 x IQ 9 (cm- 3 ) (2.4) 

K(cm-3)]J 

Eqn.(2.3) is vahd only if > 1 and 0.05 < A^ric/lO 16 ) 7 < 0.5. Using our measured 
temperature of 5.49 eV and density ~ 10 17 /cm 3 (at 4 mm), we get the number of electrons 
in the Debye sphere equal to 70 and the parameter Aj(n f /10 16 ) 1 equal to 0.06. Thus the use 
of Eqn.(2.3) for density measurement is justified. The ionic contribution to the broadening 
is observed to be much smaller and can be neglected in a typical laser ablated plasma 36. 
187], Eqn.(2.3) simplifies to 

AA = 2K (Wi) (A) 

We have used Eqn.(2.5) for estimating electron density of C II transition (3P 2 P° - 4s 2 S) 
at 392.0 am. Electron density using 1.06 jum laser wavelength and irradiance of 1.22 x 10 n 
Watt/cm 2 is found to be 10 16 -10 17 cm -3 . 

2.1.3 Plasma Expansion Velocity 

To calculate the expansion velocity of the plasma species, the temporal profiles of various 
transitions of interest were recorded at different distances from the target surface. The 
velocity of the plasma species is calculated from the delay of the peak of the signal at 
different distances with respect to the trigger pulse [39]. The velocity of CHI species at 
distance of 4 mm from the target surface in vacuum is estimated to be 1.76 xlO 6 cm/sec. 


2.1.4 Vibrational Temperature 


The vibrational temperature was estimated using C 2 Swan band head intensities [188]. 
the relative population in each vibrational level N u / can be obtained using a theoretical 
Franck-Condon factor [189] F u v> and the band head intensity 1“^" observed at frequency 


N u , = 


Cl^t 

em 


D u F u f u "V 4 


( 2 . 6 ) 
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Figure 2.2: A schematic sketch of the probe. 


where C is a constant and D u represents the correction factor for the detecting system, 
u' and u" are the vibrational levels of the upper and lower electronic states respectively. 
Since in thermal equilibrium the population of the initial state is proportional to exp 
[- G(u') h c / ke j , we have 

n < r-' u " 

■ = J ^ = (Constant)exp 


-[<?(«') - GM] 


k B T 


nb 


(2.7) 


where G(u') and G(w') is term value of the upper vibrational level u' and w', T tni) is rhe 
molecular vibrational temperature. Plotting the log of the ratio of band head intensities 
of two transitions against the difference of term value of the upper vibrational level, gives 
vibrational temperature. The vibrational temperature for the C 2 molecule is estimated to 
be 1.96 eV. 


2.2 Ion Probe Diagnostics 


The ion probes are simple and extremely useful for in-situ diagnostics of the flowing plas- 
mas. Ion probes provide the local information on plasma conditions unlike emission or 
absorption spectroscopy, laser-induced fluroscence or imaging techniques which are useful 
along the line of sight [100]. The schematic of the experimental arrangement for the ion 
probe diagnostics is shown in Fig. 2.1, where the port which is at an angle of 45° of the 
interaction chamber with respect to laser entrance port was used for ion-probe diagnostics. 
Langmuir probe consisted of a tungsten wire of 0.6 mm diameter and a power supply. Fig- 
ure 2.2 shows a schematic sketch of the probe. The probe to target distance could be varied 
continuously. The ion probe was connected through a telfon shielded wire to the power 
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supply outside the vacuum chamber. The power supply consisted of a variable -180 V to 
4-180 V DC supply. Output signal due to the electron / ion current from the probe was fed 
to a storage oscilloscope through a 50 Q terminator. The digitized ion probe signals were 
fed to the PC for further data processing. The peak of the ion probe signal (I p ) was plotted 
against the probe voltage (V p ) to get the Ip-V p characteristics. The electron temperature 
and ion density were calculated using the I p -V p characteristics of the ion probe. 


2.2.1 Electron Temperature 


The evaluation of plasma parameter from the ion probe characteristics depends on the 
shape of the probe, probe dimension and collision length etc. Ass uming Maxwellian ve- 
locity distribution of the charge carriers, the total probe current (Ip) flowing through the 
probe is given by [190]. 

/ P v \ 

( 2 . 8 ) 


Ip = /_ 0 + Ieoexp 


ksTj 


Where I +0 , fro, V p , e, T e , represent the ion saturation current , electron saturation current, 
applied probe voltage, charge on the electron and electron temperature respectively. Below 


electron saturation current, Eqn.2.8 can be used to estimate T e by using the formula, 

-dln(Ip-I +o y~ l 


T e = 


dV v 


(2.9) 


2.2.2 Ion Density 


The knowledge of the variation of the ion current with applied probe voltage offers the 
possibility to estimate the plasma ion density. Several ways of estimating ion density 
using ion probe characteristics are available in the literature [190]. The derivative of the 
square of the probe current with respect to the probe potential yields a straight line for 


the cylindrical probes, the slope of the straight line gives the density, 


Ni = 


A-xm e {-dI p 2 /dV p ) 


, i 


ZA%e 2 


( 2 . 10 ) 


where N* represents the ion density, Ap is the area of the ion probe, e and nq are the charge 
and mass of the electron respectively. 
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2.3 Fast Photography 

Fast side on views of the expanding plasma plume were made by recording the overall 
visible emission from the plasma plume with a gated intensified CCD camera system 
(ICCD I76G/2, Princeton Instruments Inc.). The complete system comprises an ICCD 
detector. ST-130 controller, PG 200 pulse generator. The controller is interfaced directly 
to a Personal Computer. 

ICCD detector : A 50 mm Nikon telephoto lens attached at the detector front is used 
to image the region of plasma. The imaging field of view and the different ports of ICCD 
detector are shown in Fig.2.3 (a) and (b) respectively. The detector has the following 
distinct sections. The front enclosure contains the image intensifier and the CCD. The 
CCD is seated on a cold finger that is in turn seated on a thermoelectric (Peltier effect) 
cooler, driven by closed loop proportional control circuitry. The back enclosure contains 
the heat removal block where coolant is circulated. A thermal sensing diode attached 
to the cooling block of the detector monitors its temperature. Heat generated at the 
exhaust plate of the cooler is conducted to the enclosure of the detector where it is carried 
away by water or other coolant. ICCD is never run without water or circulating coolant. 
Continuous flushing of the detector with nitrogen or dry air is required for ICCD while 
they are cooled. Two copper and brass fittings are used to connect the coolant ports on 
the sides of the detector. Either port may be used as the inlet point. Detector is not 
operated if the relative humidity in the lab is greater than 40%. 

The detector consists of a CCD device optically interfaced to the microchannel plate 
image intensifier through an optical-fiber coupler. This fiber optic coupling is very sensi- 
tive. With a fiber coupled detector the combined high gain of the image intensifier and 
the low read out of the CCD array results in a detector capable of responding to a single 
photoelectron. The micro chann el plate (MCP) intensifier prevents photoelectrons from 
reaching the phosphor / CCD array until a gate voltage is applied. In the shutter mode 
for the ICCD, the photo cathode potential is held approximately 200 volts more positive 
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and therefore the photoelectron image can not reach the MCP. At an exact predetermined 
time (with respect to triggering onset) the pulse generator produces a negative pulse of 
approximately 200 volts amplitude. This pulse when applied to the photocathode renders 
it (for the duration of pulse) more negative than MCP and therefore the "photo electron 
image" is transmitted to the MCP and intensified. MCP is composed of more than 1(> 6 in- 
dividual miniature electron multipliers with an excellent input to output spatial geometric 
accuracy. Intensified gain is varied by adjusting the voltage across the MCP or the voltage 
across the MCP output and the phosphor. The second parameter is factory adjusted as it 
affects both the gain and the resolution of the intensifier. 

To get an image from the CCD detector to appear on the computer screen, it must 
first be read, digitized and transferred to the computer. A block diagram of the path of the 
image signal is shown in Fig.2.4. CCD is viewed as a two dimensional grid of individual 
photo diodes (pixels, 576 X 384) each connected to its own storage bucket (called wall). 
Each pixel senses the intensity of light falling on its collection area and stores corresponding 
amount of charge in its associative bucket. These buckets are shifted to the bucket of 
shift register and then to measurement circuitry. ICCD’s can get destroyed if exposed to 
excessive light levels, higher than twice the ADC saturation levels, so care is taken not to 
expose it to room light. 

PG-200 programmable gate pulse generator : The PG-200 programmable gate 
pulse generator combines the functions of a high voltage gate pulse generator and a sophis- 
ticated timing source / delay generator. It is primarily intended to provide high speeding 
gating of proximity- focused MCP image intensifier. It is interfaced to the ICCD series 
detector equipped with internal 1 Mfi terminator. This instrument operates at high volt- 
ages. 

ST-130 controller : The interfacing of the ICCD detector and computer is through 
the controller. It controls all the timing operations. During read out, an analog signal 
representing the charge of each pixel is sent to the controller and digitized as shown in 
Fig.2.4. The number of bits per pixel is based on both the hardware settings and the 
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Computer 


Figure 2.4: A block diagram of the light path in the detector system. 
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Figure 2.5: Typical ICCD 
pulse. 

settings programmed into the controller. Recorded images are analyzed using the supplied 
software. 

Plume expansion was recorded utilizing the overall visible emission from the plasma 
plume with a gated, intensified CCD camera system. Light was collected from throughout 
the plasma volume to form a weighted two dimensional image of the pl um p, intensity. As 
explained above, MCP prevents photoelectrons from the photo emissive detector surface 
from reaching CCD array until a gate voltage pulse is applied. Fast gate pulses (acting 
as the shutter for the camera) were supplied at variable times after the firing of the laser 
for "Stop action" images of the plume expansion. Gate width was fixed at 20 ns implying 
that the electrons from the detector can reach the CCD array for 20 ns. Delay time is 
the time after the firing of laser pulse at which the gate voltage is applied. An important 
experimental parameter is the exposure time chosen for the ICCD ima ges he. the time 
for which the CCD is exposed to the emitted radiation. Over exposure results in loss of 
details in bright regions of interest while under exposure misses the fainter regions. 

Measurement of the front edge of the laser plume as a function of time helps one 
estimate the initial velocity and the stopping distances of the plume [104]. The software 
provided with the camera was tailored to calculate maximum, minimum , total and mean 
intensity of the plasma plume. The colour code helps observe different intensity regions of 
the plume. Figure 2.5 gives a typical ICCD image of carbon plasma after 200 nsec from 
the time of ablating pulse. 
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Chapter 3 


Numerical Modeling of Laser Matter 
Interaction in Presence of an 
Ambient Gas 


3.1 Introduction 

The study of laser-matter interaction has been a subject of paramount importance in 
recent years. There have been reports on laser-matter interaction in the presence of an 
ambient gas [60, 92], electric field [191, 192] and magnetic field [193]. Various analytical 
models for the plasma expansion in ambient gas have been developed [44] but they fail to 
predict details about the plasma behavior because they do not take care of the many impor- 
tant processes (for example Bremsstrahlung radiation, heat transport due to suprathermal 
electrons etc.) and hence a detailed theoretical description is necessary. The genesis of this 
study is to understand the interaction of a high power pulsed laser beam with a solid target 
in presence of an ambient gas and to develop a comprehensive numerical model for the 
°The work presented in this Chapter is based on : 

1. "Numerical modeling of laser matter interaction in the presence of an ambient gas." A. Neogi and R. 

K. Thareja, Pramana 50, 63 (1998). 
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process. Computer simulation studies help optimise the plasma parameters essential for 
various technological applications and studying astrophysicai processes without actually 
doing the experiment. 

The present Chapter deals with the case of intermediate intensity regime. 10 h -10 r - 
Watt /cm 2 of laser power and high ambient gas pressure. We have used a modified version 
of the laser fusion code MEDUSA [51], MED101. This is an one di m ensional Lagrangian 
code which calculates hydrodynamic and thermodynamic parameters of the laser produced 
plasma in vacuum. MED 101 is a well developed code which can take care of the double 
and triple layered targets and can calculate X-ray laser gain. We. have made modifications 
[58] in this code to implement it for ambient gas. The Rankine-Hugoniot relations for 
shock waves are used to calculate the changing boundary conditions that is temperature 
and pressure. The Navier Stokes equation is changed due to Reynolds stress arising due 
to turbulent boundary layer mixing and at later stage due to diffusion and interparticie 
collisions. The thermodynamics of the system is also changed accordingly. The equations 
are solved using finite difference method. 

The modified code MED 101 is used to simulate the experimental results where the 
laser produced carbon plasma expanded into a background of He, Ne, Ar and Xe gas. 

3.2 Numerical Model 

3.2.1 In Vacuum 

A brief outline of the one dimensional Lagrangian code MEDUSA is presented below. A 
detailed description and the mathematics of the original code is given by Christiansen et 
a 1 [51]. The plasma is assumed to consist of a charge neutral mixture of electrons and 
various species of ions. MEDUSA treats the thermodynamical properties of the electrons 
as one subsystem with internal energy U e per unit total mass, temperature T e , pressure 
P e and specific heat ratio -y e and the ions as second subsystem with a corresponding set of 
variables. Electric and magnetic fields are neglected and the two subsystems are coupled 
together by a common velocity and by energy exchange due to collisions between them. 
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This is the one fluid two temperature model, the whole of the plasma is considered as one 
fluid in which electrons and ions have different temperatures. The internal energy per unit 
mass is denoted by U (subscript e and i are omitted), 


U = 


PV 

(7 - 1) 


(3.1) 


where P is pressure, V, the volume per unit mass and 7 the specific heat ratio. The two 
subsystems have independent equations of state 

U = U(p.T). P = P(p. T) (3.2) 


where p is the density of the plasma. 
According to the first law of thermodynamics. 


dQ_dU_ dW 
dt dt dt 


where Q, U and W represent heat, internal energy and work done. 




dT , dp dV 

l . / 3 r — + p — 

dt 1 T dt dt 


(3.3) 


(3.4) 


Cv ~ { 6T } ; 


0T = (- 7 —) 
op T 


The source terms Si and S e for ions and electrons respectively are written as; 

Si = - R + Q s 


(3.5) 


(3.6) 


S e = H e + R + L+X (3-7) 

where H represents the flow of heat due to thermal conduction, R is the rate of energy 
exchange between ions and electrons, L the rate of Bremstrahlung emission. X the rate 
of absorption of laser light and Qs is the rate of viscous shock heating. In plasma, if 
shock wave is formed, thermal conduction or energy transport is due to electrons whereas 
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viscosity or momentum transport is due to the ions. The ions are heated up due to 
viscous pressure of the shock wave which is represented by Q s . the viscous shock heating. 
Mathematical expression of each term appearing in Eqn.(3.6) and (3.7) are given below 
(in SI unit) 

(1). The flow of heat due to thermal conduction ( H ) in Eqn.(3.6) and (3.7) is given by 

H = -V.kVT (3.8) 

P 

where k is thermal conductivity. Expressions for thermal conductivity due to ions (*,) and 
electrons ( K e ) are given by [193] 

= 4.3 x IQ" n T? {log A) " 1 M"5 2 - 2 (F)~ 1 (3.9) 


= 1.83 x lO~ m T}{logAy l z{z 2 ) 


(3.10) 


where Ti and T f , are the ion and electron temperature, log A is the Coulomb logarithm 
given by Eqn.(3.11), M and z are the average mass and charge number given by Eqn.(3.12) 
and (3.13) respectively 

_ 1 

A = 1.24 x 10 6 Tj— (3.11) 

z 

M = Y, ( 3 - 12 ) 


z = YL fkZk 


(3.13) 


where M fc , Zk are the mass and charge number df the individual species. 

(2). The exc han ge of energy between ions and electrons ( R ) in Eqn.(3.6) and (3.7) 

depends on the difference of temperature between them and is given by 

R = -u ex n e k B (T i -T e ) (3.14) 

P 

where u ex is energy exchange rate given by [194] ^ 

M~ l z 2 e i n i logAml 3 

= - 7=— 0 ( fc S T e ) - 

32v27re5mff 
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where m P and m# represent electron mass and proton mass respectively. Substituting for 
Ufx in Eqn.(3.14) we get 

R = 0.59 x 10 _8 n e (Ti - T e )TpJF^logA (3.16) 

(3) . The plasma loses heat mainly through Bremsstrahlung radiation. This ocuurs when 
electrons come close to the nucleus, get accelerated and emit photon. The Bremsstrahlung 
radiation emitted by the electrons is represented by L in Eqn.(3.7). The radiation is 
assumed to escape from the plasma and is thus an energy loss (L < 0) term. For a 
Maxwellian gas L is [194], 

T = - w (3.17) 

2A7r€o 3 c i m e m f fMh 

Inserting for v e , the electron velocity, Eqn.(3.17) becomes 

L = -8.5 x \Q- u n e T0M~ x (3.18) 

(4) . Absorption is assumed to occur via inverse Bremsstrahlung at densities below the 
critical density. If P i (Ro) is the laser power incident on the plasma boundary at r = Ro 
then at any point r inside the plasma from the boundary we have 

P L (r,t) = exp[-n(i?o - (3.19) 

Here a is the absorption coefficient defined earlier [Eqn.(1.2)j. The term X(r, t) in Eqn. (3.7) 
is defined as laser power absorbed per unit mass, 


X(r,t) 


T’z.Cr.t) 

dM 


(3.20) 


where d M is an infinitesimal mass element. 

(5) . The ions are heated up due to viscous heating caused by the shock wave, the viscous 
heat term is given by [19] 


Qs 



( 3 . 21 ) 
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where q is a viscous pressure given by 


dv 

dt 


< O(Compression) 


(3.22) 


5 = 0, 


av_ 

dt 


> O(Expansion) 


(3.23) 


where b is a coefficient of viscosity usually chosen as b = 2 1/f2 . 

(6). The plasma is considered to be a fluid, the motion of the plasma is governed bv the 
Navier-Stokes equation, 


dv 


P 


dt 


= -VP 


(3.24) 


where P is the total pressure and v is the velocity given by 

. , dr 


(3.25) 


(7) It is necessary to specify both hydrodynamic and thermodynamic boundary conditions 
at the inner and outer boundaries. For the inner boundary we choose 

. v (r — 0) = 0. 


At a movi ng outer boundary point r=Ro we distinguish the following four cases. 

1. If pressure at the external boundary is zero and no thermal flux, usually referred to 
as zero thermal flux case, 

P (Bo) = 0 

2. The velocity of the external boundary is zero, ion and electron temperature of exter- 
nal boundary is changing as some function of time, 

v (Bo) = 0, T, (R c ) = T t (t), T e (R 0 ) = T e (t). 

3. The pressure applied at the outer most boundary is some function of time and no 
thermal flux (zero thermal flux), 

P (Re) = P(t) 
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4. The velocity of the external boundary is changing as some function of time and no 
thermal flux (zero thermal flux), 
v (Rq) = v(t) 

Case 1 is used for general problems involving the absorption of laser light, case 2 is used 
to study simplified thermal conductive problems, and cases 3 and 4 are used to study 
simplified hydrodynamics problems e. g. compression and expansion of a plasma. 

These equations are solved using finite difference method. The whole length of the 
plasma is divided into finite number of meshes and time is divided into large number 
of small time steps. The plasma parameters like position of the plasma front, velocity, 
temperature, pressure and density etc. are calculated as a function of time and space 1 - The 
original code which calculates the plasma parameters in vacuum is modified for the case 
of plasma expansion in high pressure ambient gas and is described in detail in the next 
section. 

3.2.2 In Ambient Gas 

The plasma expanding in ambient gas can broadly be divided into two regions in accor- 
dance with external gas pressure. (1) If the external gas pressure is low such that the 
mean free path of the plasma particles exceeds or is comparable with the dimension of 
experimental chamber (~ 10 cm). In this case shock wave does not form and there is a 
strong interpenetration of laser plasma and the ambient gas [60]. The calculated shock 
thickness is less than the mean free path of the plasma particles [98]. (2) The external 
pressure is high (say > 1 Torr) and the mean free path is much less than the dimension of 
the exper imen tal chamber. In this case shock wave forms in the ambient gas, the calculated 
shock thickness is considerably greater than the mean free path of the plasma particles. 
The present Chapter deals with the physics of plasma expansion in high pressure external 
gas and the formation of shock wave. 

When the velocity of fluid in motion (or any disturbance) either becomes comparable 
with or exceeds that of sound (supersonic fluid) it gives rise to shock wave in the ambient 
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Figure 3.1: (a) The moving shock front as seen from reference frame fixed to earth; (b) 
The stationary shock front as seen from reference frame moving with the shock front. 

medium. Mathemetically, we can define a shock wave as any abrupt disturbance that prop- 
agates through the gas causing a change of state. The laser produced plasma expanding 
in high pressure ambient gas giving rise to shock wave in the medium is a typical example 
[92], 

Figure 3.1(a) shows the propagation of a plane shock wave propagating from left ro 
right in a medium with speed vj + v s assuming the fluid ahead of and behind the shock 
wave in a steady state described by the state variables in the figure. The flow velocity 
ahead of the shock is taken to be v !; and the gas behind the shock is then set into motion 
with a velocity v 2 . If Vi is set equal to zero it becomes a special case where flow velocity 
ahead of the shock is zero and shock moves with velocity v s . We will solve for general case 
when v x ^ 0. Our goal is to determine the gas properties behind the shock wave such as 
density p 2 and the flow velocity v 2 . 

To analyze the shock wave it is convenient to shift to a coordinate system moving 
along with the shock as shown in Fig.3.1(b). The conservation equations for any control 


volume are 
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The continuity equation : 



( 3.26) 

The momentum equation : 

|(pt.)+£(p+p» 2 )=°. 

(3.27) 

The energy equation : 

9 pv 2 d [ P 1 

S (pl + T ) + Tx r (E + 7 ) j = °- 

(3.28) 

Where , U = C V T 

(3.29) 

and H = U + — 

P 

(3.30) 

Equation of state (for perfect gas) : 


P = pBT. 

(3.31) 

For steady state flow the time derivative vanishes and the Eqns.(3.26 - 

3.28) reduce to 

Pi«! = P2*4 

(3.32) 

Pi + pWi = Pi + P 2 V 2 

(3.33) 

2 ui 2 

Ei + — = H 2 + — 

2 2 

(3.34) 

P = pRT 

(3.35) 


where v^ and v^ are the velocities in the moving frame attached to the moving shock. In 
the stationary coordinate system flow velocities vi and v 2 are different from V T and v, . 
The thermody nami c properties such as pressure, density, temperature and speed of sound 
are not influenced by the transformation of coordinate system. The relations between the 
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flow velocities in the two transformations can be expressed as 

v[ = v , 


(3.36) 


Vo = V, + (|>1 — Vo). 


(3.37) 


The flow Mach number, Mj is defined by 


— = — = M„ 

a\ a\ 


where aj represents the velocity of sound. From the relation Eqns.(3.32 - 3.35) and 


Eqns.(3.36 - 3.38) and using the relation a! = 7 RTinP = pRT one obtains 


vo - v x ,2 1 

a, 7 + 1 V M s 


(3.39) 


(7 + 1)' 


(7-1) 
(7 + 1) 


(3.40) 


P 2 _ (7 + 1 )M„ 2 

Pi “ (7 — 1)M S 2 + 2 


(3.41) 


T 2 r2 7 Mj 2 — (7 — 1)1 r. 3 1 

■ Ji-l'h+prj l (7 “ 1) "' +2 ‘ 

For strong shocks, M s >> 1 , Eqn.(3.41) reduces to 

P 2 ^ 7 + 1 

pi 7-1 


(3.42) 


(3.43) 


As the plasma moves at a velocity of 10 6 - 10 7 cm/sec into an ambient gas it acts as 
a supersonic piston giving rise to a shock wave which moves at a faster velocity than the 
piston [16]. For high pressure gas there is a constant discontinuity between the piston and 
the swept up gas and there hardly is any intermixing of the two. A rough sketch of the 
process is given in Fig.3.2. Region I, the undisturbed gas where pressure and temperature 
are P 0 and T 0 , region II is the compressed gas ahead of the plasma front and the region 
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Figure 3.2: Region III is expanding plasma, region II is compressed gas and region I is the 
undisturbed gas shown in the reference frame fixed to the earth. P 0 / (P c ) and T 0 / (T c ) 
refer to pressure and temperature in region I and II respectively. 

HI is the plasma region. Between region III and II is the plasma front which is referred 
to as supersonic piston. Between region I and II is the shock front which is moving faster 
than the plasma front. 

The velocity of the plasma front v v and the velocity of shock front v, are related by 
the expression [16] 


(7 + 1) 


v r + 


'(t+JF 

16 


Up 2 + 0L\ 2 


(3.44) 


The boundary condition on the region IH is defined by the pressure and temperature 
of the region II. Here the pressure (P c ) and temperature (T c ) in compressed gas region II 
across a stationary shock front are given by (from Eqs. (3.40) and (3.42)) 


Pc = Po 


2 7 


(7 + 1) 


M? — 


(7-1) 


(7 + l)j 


(3.45) 


T C =T 0 


27M, 2 — (7 — 1) 


(7 + 1 w 


[fry — DM, 2 + 2 ] (3.46) 
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It follows from Eqn.(3.45) and Eqn.(3.46) that the boundary condition is changing 

with time, the changed boundary conditions are 
P(Ro) = Po if M, < 1 
= P c if M, > 1 
T e (Ro) = To if M, < 1 
= T c if M, > 1 


(3.47) 


If the Mach number is less than or equal to unity then shock waves are not formed 
and the boundary conditions are defined by Po and To but if Mach number is greater than 
unity then shock waves are formed and pressure and temperature are given by P c and T c 
respectively. 

The boundary condition for the temperature is applicable to the electron temperature 
only. It is because of the higher concentration of electrons at the plasma front due to their 
higher mobility compared to ions. The velocity is being calculated at every step. That is, 
after every small time step the code calculates the plasma front velocity and simultaneously 
calculates the shock front velocity and the corresponding pressure at the boundary and 
with that pressure it calculates the next plasma front velocity and the process continues. 

When plasma moves into an ambient gas it is almost similar to the case of hot jet 
coming out into ambient atmosphere. The moving plasma undergoes turbulent boundary 
layer mixing due to which it experiences viscous force on its boundary. According to 
Prandtl’s mixing length theory, the total Reynolds stress due to turbulent mixing in the 
boundary layer is given by the expression [195]. 


PO^r c — p0% 1 b\ [Umax t^min] £ 

by by 


(3.48) 


where po is the density of ambient gas, e r is kinematic viscosity, iw is the velocity of the 
center of flow and v m i„ is at the boundary, xi is an empirical constant to be determined 
from experiment, bi is the width of the m i xin g zone. 

At the outermost boundary, velocity is zero, so v m { n = 0 and writing tWr = v, Eqn.(3.48) 
simplifies to 
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r = pQX\V~ 

The shear force experienced by a particular mesh will be 


(3.49) 


Fdrag = t(61) (3.50) 

51 is the length of a particular mesh which is varying with time. The Navier-Stokes equation 
( Eqn.(3.24) ) changes to 


p— = -Vp - t(AZ) (3.51) 

As the plasma front expands into the ambient gas the pressure just inside the boundary 
decreases and the pressure outside increases. At a later time the pressure inside the plasma 
front is comparable to the pressure outside. At that time the piston becomes leaky and 
outside gas diffuses into the plasma giving rise to collisions due to which the expanding 
plasma experiences a large resistive force, 

Pint « 

or 

Fdiff — c 2 (^—) Po (3.52) 

Pint 

p ext is the pressure just outside the plasma front and pi nt is the pressure just inside the 
plasma front. On addition of a diffusive term the modified Navier-Stokes equation reads 

p^ = -VP - r(AZ) - F^. = -VP - pox.v 2 - c 2 (§^) 1 po (3-53) 

at -Lira 

This equation contains three constants viz. aq, c, and c 2 to be determined from the 
experiments. 

The thermodynamics of the system changes in the presence of the ambient gas. The 
kinetic energy (KE) lost due to drag and diffusion is expressed as heat which increases the 
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temperature of the expanding plasma. Some heat is also lost to the gas . For the sake of 
simplicity it is assumed that the heat lost to the surrounding gases is negligible, loss in 
KE is expressed as heat is. 

^m[u, 2 -v f 2 ]= mC p 8(J c ) (3.54) 

J c is joules constant, C p is specific heat at constant pressure and 9 is the temperature 
difference. The heat lost to the shock wave and the work done in compressing the gas 
has been taken care of by the boundary conditions because it changes the pressure and 
temperature inside the plasma. 

3.3 Results 

With the modifications as above viz. including the viscous and diffusive force in the Navier 
Stokes equation and changing the boundary conditions, the simulation is carried out (with 
an initial step of t=10~ 18 sec having mesh size of 60) and matched with the experimental 
data [93]. Here graphite target was ablated using KrF laser (A = 248 am, 30-ns FWHM, 
2x2 mm 2 spot size on the target, fiuence=20 I /cm 2 ). The ablation was carried out in 
four different buffer gases at a pressure of 300 Torr. The temporal evolution of the plasma 
front as it propagates in He, Ne, Ar and Xe gas is shown in Fig.3.3(a) - 3.3(d). with a 
dotted line. One of the simulated curves for a particular gas (say Xe) is compared to the 
experimental curve and constants xi, Cj and Co in Eqn.(3.53) are determined. It was found 
that xi = 2.0 x 10 3 /m 2 , c\ = 1 and c 2 = 1.0 x 10 9 /m 2 . These constants once determined 
are then used to get other simulated curves (in case of He, Ne and Ar) which are compared 
with respective experimental curves. The simulated results shown in solid line match 
well with the experimental results, shown in dashed line till 300 nsec. The plasma front 
propagates faster in lighter gas (like He) than denser gas (like Xe). It is interesting to note 
that after 300 nsec the simulated value is more than that of the experimental one. This 
is because after 300 ns the diffusion increases to such an extent that a reflected shock is 
generated within the plasma plume which moves in the opposite direction, towards the 
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Figure 3.3: The position of the plasma front versus time in four different ambient gases 
(a) He (b) Ne (c) Ar and (d) Xe. Experiment ( ) Simulation ( ). 
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Figure 3.4: The simulated temporal profile of velocity of the plasma front in Xe gas. The 
inset shows the oscillations observed at around 300 ns on an enlarged scale. 

target surface, dam ping the motion of the plasma as observed in the experiment. If there 
is mixture of fluids whose composition is different at different points then the equations 
of fluid dy nami cs have to be considerably modified [196]. A diffrent set of equations are 
needed to model this part which are not considered by us. 

Figure 3.4 shows simulated temporal profile of plasma front velocity with Xe gas as 
the background. The velocity decreases with time but at around 290 ns fluctuations in 
velocity are observed. Figure 3.5 shows the intensity of the shock wave going into the 
plasma as a function of t im e as observed in simulation. The shock intensity is maximum 
at around 290 ns, that is at the same time when velocity suddenly drops and increases 
again. Experimentally it was observed that a shock wave had gone into the plasma front 
at around 300 ns causing instability in the velocity of the plasma front. The simulated 
results match very well with the experiment. It has also been observed in simulation that 
the intensity of shock wave going into the plasma in case of He is orders of magnitude 
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Figure 3.5: The simulated temporal profile of shock intensity at the plasma front in case 
of plasma expanding in Xe gets. 

lower than that of Xe. It is because Xe being a denser gas gives a greater thrust to the 
moving plasma compared to that of He. Figure 3.3(d), that is position of the plasma front 
expanding in Xe and Fig.3.4 that is the velocity of the plasma front in the same ambient 
condition is replotted in Fig.3.6, with a common time axis, position and velocity axes are 
shown along the Y axis in left and right respectively. The whole graph is divided into seven 
different regions to show that slope of the position graph obtained by simulation gives the 
simulated velocity graph. For e.g. in region (f) the slope of the position curve suddenly 
tends to 0, the velocity also tends to 0 at that point as seen in the velocity curve whereas 
in region (d) and (g), the velocity is constant. 

The simulated temporal profile of temperature, density, pressure and electron density 
of the plasma front in case of Xe (dashed fine) and He (solid fine) as buffer gas are shown 
in Fig.3.7 - 3.10. Xe being denser than He gives higher resistance to the moving plasma 
resulting in higher temperature, pressure, density and electron density of the plasma front. 
The temperature of the plasma front attains a maximum value at 30 nsec (Fig.3.7) which is 
the pulse duration of the laser used and then decreases due to the cooling effect of the inert 
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Figure 3.6: Fig.3.3(d) and Fig.3.4 are replotted with a common time axis. 
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Figure 3.8: Density of the plasma front expanding in Xe and He as obtained from simu- 
lation. Dotted line is for the case of Xe and solid line for He as ambient gas. The inset 
shows density of the plasma front expanding in He on enlarged scale neglecting the initial 
points. 
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Figure 3.9: Pressure of the plasma front expanding in Xe and He as obtained from simu- 
lation. Dotted line is for the case of Xe and solid line for He as ambient gas. 


gases. The density profile in case of Xe (Fig.3.8) shows a sudden rise at around 290 nsec 
due to shock wave coming into the plasma . However, rise in density at around 290 nsec 
(Fig.3.8) and fall in its temperature at the same time (Fig. 3.7) give no significant change in 
pressure as expected from the perfect gas equation of state. Figure 3.9 shows a simulated 
pressure profile at 290 ns. The inset in Fig.3.8 and Fig.3.10 show the density profiles 
and electron density profiles for He on a magnified scale. It shows that electron density 
increases similar to that of Xe but increase is gradual and low in case of He compared to 
Xe where change is abrupt and high because of its higher density. 

Figure 3.11 shows the effect of different terms which are added in the Navier-Stokes 
equation other than the boundary conditions. The viscous drag term plays the major role 
in decelerating the plasma front. The diffusive term comes into play at around 300 ns and 
is much less compared to viscous drag. 

We applied our model to another experiment where excimer laser produced copper 
plasma was expanding in Ar, He, N 2 and H 2 at 5 mbar ambient pressure. Figure 3.12 gives 
the simulated temporal profile of electron density for four different gases. It is observed 
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Figure 3.10: Electron density of the plasma front expanding in Xe and He as obtained 
from s imu lation. Dotted line is for the case of Xe and solid line for He as ambient gas. 
The inset shows electron density of the plasma front expanding in He on enlarged scale 
neglecting the initial points. 



Figure 3.11: The simulated position of the plasma front as obtained from (a) Eqn.(3.24) 
i.e. no change in Navier-Stokes equation, (b) Eqn.(3.51) i.e. on adding viscous drag term 
only, (c) Eqn.(3.53) i.e. on adding viscous drag term and diffusive force term. 
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Figure 3.12: The electron density of the plasma front expanding in Ar, N 2 , He, H 2 as 
obtained from simulation. 

that electron density is highest in the case of Ar followed by N 2 , He, and H 2 in agreement 
with the experimental observations [93]. 

Finally we applied our model to our own experiment in which laser ablated carbon 
plasma expanded in 100 mTorr air where the shock wave just forms. It is found that the 
position of the pl asm a front obtained from IOCD photographs matches very well with the 
simulated value for 100 mTorr air as shown in Fig. 3.13. The temperature of the plasma 
plume in vacu um matches well with the simulated value as shown in Fig. 3.14. 
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Figure 3.13: The position of the plasma front versus time in 100 mTorr ambienr air. 
Experiment ( ), Simulation ( ). 



Figure 3.14: The temperature of the plasma plume versus 
), Simulation ( -)• 


time in vacuum. Experiment (- 



Chapter 4 


Laser Produced Carbon Plasma in 
Low Pressure Ambient Atmosphere 


4.1 Introduction 

A high temperature and high density transient plasma produced by focusing the out- 
put from a Q-switched laser onto a solid target has attracted a great scientific concern 
due to its recent technological applications, in particular to material research [10], fluid 
dynamics [196], cosmology [197], chemical physics [112] and many other fields. Laser ab- 
lated plasmas of metals in presence of background gas have been used in the recombination 
phase for the generation of laser oscillations [198, 199] and as strong x-ray and vacuum 
ultra-violet sources [200, 201]. The interaction between the laser plasma and ambient 
gas provides information on collisional, collective and electromagnetic processes in astro- 
physical and laboratory plasmas [8]. Since the background gas effectively controls the 
plasma par am eters important for numerous technological applications, the knowledge of 
°The work presented in this Chapter is based on : 

1. "Dynamics of laser produced carbon plasma expanding in low pressure ambient atmosphere" , A. 

Neogi, A. Mlshra and R. K. Thareja, J. AppL Phys. 83, 2831 (1998). 
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the dynamics of the plasma plume expanding in ambient atmosphere is indispensable. 

There have been several reports regarding the propagation of laser produced plasmas 
in ambient atmosphere, however, the dynamics of carbon plasma expansion in ambient gas 
at low pressure is still not clear. In the present Chapter a detailed study of laser produced 
carbon plasma expanding in vacuum (0.1 mTorr) and low pressure ambient atmosphere (5 
mTorr and 100 mTorr) is presented using ICCD photography and time resolved emission 
spectroscopy. The plasma plume is observed to break into fast and slow components in 5 
mTorr air. Double peak structure is observed in the temporal profile of Cl, CII and CIII 
transitions. Probable reasons for the stratification of the plasma and double peak behavior 
are discussed. 

4.2 Experimental Details 

The experimental set up used in the present study is described in Chapter 2. The energy 
of the laser pulse was 800 mJ/pulse, and was focussed on to the rotating graphite target 
resulting in irradiance of ~ 1.22 x 10 11 Watt/cm 2 . The plasma expanding in 0.1 mTorr. 5 
mTorr and 100 mTorr air was studied using ICCD photography and emission spectroscopy. 

4.3 Results and Discussion 

4.3.1 ICCD Photography 

Figure 4.1 shows ICCD photographs of the plasma plume at three pressures 0.1 mTorr, 
5 mTorr and 100 mTorr at various delays with respect to ablating laser pulse with 20 
ns exposure time. Each image is normalized to its maximum intensity. In vacuum the 
images were visible till 500 ns and became highly diffused beyond 500 ns. On increasing 
the pressure to 5 mTorr a significant increase in intensity was observed. At 600 ns the 
plume breaks into two components, fast and slow. The velocity of the two components 
is es tima ted to be 1.0 xlO 6 cm/sec and 7.5 x 10 5 cm/sec and hence termed as fast and 
slow components. The fast component Is bulky consisting of the major part of the plume 
and the slow component is the minor component as shown in Fig. 4.1. This trend was seen 
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in all other subsequent photographs. However, at 1500 ns the slow component became 
bulky with the fast component becoming the minor component. At around 2000 ns only 
one component is observed. At later times the plume remains as a single component 
of an irregular shape similar to that observed by Anan’in et ai [25] . At 100 mTorr the 
stratification of the plasma plume is not observed. The dim ensions of the pl um e get 
reduced due to pressure of the ambient gas from outside. 

4.3.2 Emission Spectroscopy 

The time resolved emission spectra of carbon plasma were recorded at different distances 
away and parallel to the target surface by moving the monochromator in a plane perpen- 
dicular to the target surface. The temporal profile of Cl transition 2p 3 3 D° - 8f F(5/2) at 
399.7 nm, CII transition 3d 2 D - 4f 2 F° at 426.7 nm, CIII transition 3s 3 S - 3p 3 P° at 465.0 
nm and CIV transition 3s 2 S - 3p 2 P° at 580.1 nm were recorded at different distances 
from the target surface at various pressures of air. Each of these transitions is one of the 
strongest transitions in their respective ionic states. The e mis sion lines are identified using 
the information available in the literature [184]. Velocity, intensity and nature of the tem- 
poral profile of all different species are studied in vacuum (0.1 mTorr) and low pressure 
ambient atmosphere (5 mTorr and 100 mTorr) [The Fortran program used to calculate 
these parameters is included in the appendix A]. 

The velocity of each of the species Cl, CII, CIII and CTV is calculated from the delay 
of the peak of the signal at different distances with respect to the trigger pulse [10,12]. It 
is observed that CIV ions have the highest velocity followed by CIII. Cl and CII ions in 
case of vacuum. The expansion velocity of various species at a distance of 4 mm from the 
target surface is shown in Table 4.1. 

The difference in velocity for the different ions and neutrals can be explained in terms 
of mobility which is proportional to charge to mass ratio (qj/m;). CIV ions have the highest 
mobility followed by CHI, CII and Cl. As the plasma expands in vacuum, the electrons 
strive to overtake the ions resulting in an uncompensated space charge and electric field 
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Table 4.1: Velocity of different species of carbon in 0.1 mTorr, 5 mTorr and 100 mTorr of 
air. 


Species 

Velocity (10 b cm/sec) 


P = 0.1 mTorr 

P = 5 mTorr 

P = 100 mTorr 

Cl 

1.58 

1.75 

1.62 

CII 

1.24 

1.90 

2.16 

CIII 

1.76 

1.71 

1.50 

CIV 

1.90 

1.65 

1.54 


that decelerates the electrons and accelerates the ions [77]. This results in higher velocity 
for high ionic state. The reason for the observed velocity of CII lower than that of Cl is not 
clear. It may be due to the fact that CII ions experience Coulomb repulsion on account of 
higher positive charge of CIII. The increase in pressure to 5 mTorr decreases the velocity 
of CIII and CIV ions as expected but the velocity of Cl and CII ions increases as compared 
to their vacuum value. Here CII ions have the highest velocity followed by Cl, CIII and 
CIV, as shown in Table 4.1. With a further increase in pressure to 100 mTorr, the velocity 
of Cl species decreases (still greater than that in vacuum) along with CIH and CIV ions 
but the velocity of CII ions increases. 

Figure 4.2 shows the variation in spectral line intensity of Cl, CII, CIII and CIV tran- 
sition with distance for three different pressures, 0.1 mTorr, 5 mTorr and 100 mTorr. The 
intensity of the spectral lines of the ionic species increases in presence of gas. Further- 
more, it is observed that the line intensity in case of CIII transition increases as pressure 
is increased from 0.1 mTorr to 5 mTorr but it decreases again at 100 mTorr. This decrease 
in intensity after a certain critical pressure (~ 100 mTorr) is observed for species with 
higher charge (CHI, CIV) and not for species with lower charge (Cl, CII). On the other 
hand, the intensity of CII and Cl species increases in 100 mTorr after a certain critical 
distance compared to 5 mTorr. This increase in line intensity with external gas pressure 
and subsequent decrease after a certain critical pressure has also been observed earlier and 
the reasons for the above observation have been discussed by Abhilasha and Thareja [65]. 
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Figure 4.2: Spatial variation of the emission intensity of (a) Cl, (b) CII, (c) CIII and (d) 
CIV transition in 0.1 mTorr, 5 mTorr and 100 mTorr of air. 
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Rise in emission line intensity in presence of gas (5 mTorr) is caused duo to enhanced 
collisions of the plasma species with the external gas particles. However, when pressure 
is increased (100 mTorr), the external gas acts as a heat sink which lowers the temper- 
ature and favors three body recombination. The higher species recombine to give !ower 
species which results in decrease in emission intensity for higher species (CIH and CIV) 
and increase in emission intensity for lower species (CII). 

Now we discuss the nature of the temporal profile of various transitions in case of 
vacuum (0.1 mTorr), 5 mTorr and 100 mTorr of ambient pressure. In vacuum a single 
peak structure is observed for the transition of all the species. In presence of gas ( 5 mTorr 
and 100 mTorr ) a peculiar double peak structure is observed in the temporal profile of 
Cl transition at 399.7 nm. CII transition at 426.7 nm and CIII transition at 465.0 nm 
as shown in Fig.4.3. Though the appearance of this double peak structure is reported 
due to stratification of the plasma into slow and fast components [25] but the re;isons 
for this behavior are still not clear. Moreover, many important features of this double 
peak structure are overlooked. Some of the salient features of this structure observed 
experimentally by us are summarized below. 

Double peak in the temporal profile is observed at 8 mm (from the target surface) 
for Cl, CIII transition and at 10 mm for CII transition when pressure in the chamber is 5 
mTorr. The peak which appears early in time i.e. the delay with respect to the ablating 
laser pulse is less, is termed as fast peak PI whereas that which appears later or the delay 
with respect to ablating laser pulse is more is termed as slow peak P2. No double peak 
structure is observed for CIV transition. When pressure is increased to 100 mTorr, the 
onset distance for double peak decreases by 2 mm. 

At 5 mTorr, it is found that the velocity of PI, the fast peak, is greater than that in 
vacuum where only single peak is observed as shown in Fig.4.3(a and b). The velocity of 
P2, the slow peak, is lower than the velocity of PI for the transitions of Cl, CII and CHI 
species. Figure 4.4(a-c) show the spatial variation of velocity of PI and P2 at 5 mTorr 
and 100 mTorr case for the case of Cl, CII and Cm transitions respectively. When the 
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Figure 4.3: Temporal profile of Cl transition, at 399.7 nm, CH transition at 426.7 nm, CHI 
transition at 465.0 nm and CIV transition at 580.1 nm. PI and P2 are respectively first 
and second peak (shown only for the 100 mTorr case). Dashed line ( ) is for P — 0.1 

mTorr at a distance of (a) 10 nnm (b) 10 mm (c) 12 mm (d) 8 mm. Dashed dot line (-*-• 
- • - ) is for P = 5 mTorr at a distance of (a) 10 mm(b) 10 mm (c) 12mm (d) 8 mm. Solid 

line ( -) is for P = 100 mTorr at a distance of (a) 10 mm (b)10 mm (c) 8 mm (d) 8 

mm. 






Figure 4.4: Velocity of fast (Pi) and slow peak (P 2 ) for (a) Cl, (b) CII and (c) GUI species 
at 5 mTorr and 100 mTorr air respectively, (d) velocity of CIV species at 5 mTorr and 100 
mTorr air respectively. 
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pressure of the ambient gas is increased from 5 mTorr to 100 mTorr it is found that velocity 
of PI increases and that of P2 decreases (except in CD). Figure 4.4(d) shows the spatial 
variation of velocity of CIV species. The velocity of CIV ions decreases with the increase 
in pressure of the external gas and no double peak structure is observed (Fig.4.3(d)). 

PI accelerates with distance (till 18 to 20 mm) whereas P2 decelerates for 5 mTorr 
and 100 mTorr pressure for the transitions of Cl, CII and CHI. A sudden rise in the 
acceleration of PI at 14 to 16 mm for all the three transitions is observed. 

It is observed that in case of 5 mTorr pressure, Pi is stronger than P2 till 12 mm for 
Cl, till 14 mm for CII, till 10 mm for CIH respectively [Fig.4.3(a), (b), (c)] whereas at 
distances greater than 12, 14 and 10 mm for Cl, CII and CHI respectively P2 is stronger 
than PI as shown in Fig.4.5. The figure shows the temporal profile of Cl, CII and CIH 
transition at 14 mm, 16 mm and 14 mm respectively. When pressure is increased to 100 
mTorr, PI remains stronger than P2 at all the distances for Cl and CII transition as shown 
in Fig.4.3(a) and (b) whereas for CIII transition, P2 is stronger than PI for all distances, 
Fig.4.3(c). The ICCD photographic observations are in accordance with spectroscopic 
observations where it is observed that PI is stronger than P2 at earlier distance and P2 is 
stronger than PI at later distance for all the three transitions. 

The reasons for the stratification of the plasma plume expanding in low pressure am- 
bient atmosphere as observed through ICCD photography, velocity of the different carbon 
species and double peak nature of the temporal profile of the carbon species and their 
relative intensities as observed through time resolved emission spectroscopy is explained 
below. 

As the plasma front moves in an ambient gas, at earlier distances with respect to the 
target surface, the gas in front of ic gets compressed due to its inability to move as fast as 
the plasma front. The plasma front gives a thrust on to the compressed gas which sends 
an equal and opposite compression pulse resulting in the compression of the plasma front. 
Once the compression pulse is over the thin compressed plasma front uncompresses itself 
due to which the velocity of the fast moving species closer to the plasma front decreases 
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whereas the velocity of the species in the core increases to conserve the momemtum of the 
plume. Experimentally we have observed that the velocity of Cl and CH ions increases 
whereas that of CHI and CIV ions decreases at 4 mm on changing pressure from 0.1 mTorr 
to 5 mTorr [Table 4.1]. 

A probable cause for the stratification of the plasma plume in case of 5 mTorr air mav 
be the direct collisional ionization of the air close to the plasma front by the high energetic 
plasma constituents emanating from the ablated plume. These ionized species diffuse into 
the plasma front and collide with various carbon species in the plume. Since the Coulomb 
collision cross section depends on the velocity and charge, various carbon species in the 
plume will have different collision cross section. The collision cross section (cr c ) due to 
Coulomb collision depends on the charge and velocity as [202] 


where q, the charge of carbon species and v, the velocity of the species. The species with 
higher velocity and lower charge [CII] have lower collision cross section compared to the 
species with lower velocity and higher charge [CIII, CIV]. A rough estimate of the ratio of 
the collision cross section of the species CII : CIII : CIV is 1 : 22 : 79 for the 5 mTorr case. 
Since the collision cross section determines the relative number of collision so species with 
higher collision cross section [CIII, CIV] undergo more number of collisions compared to 
[Cl, CII] (lower collision cross section). The collided species give rise to slow component 
P2 whereas uncollided species give rise to fast component PI. That is why it is observed 
that in case of CIII and CIV slow peak P2 is stronger whereas in case of Cl and CII fast 
peak PI is stronger. Due to the ongoing collision process P2 becomes stronger than PI at 
a larger distance for all the species as shown in Fig.4.5. 

In case of 100 mTorr double peak structure is observed spectroscopically but no strat- 
ification is observed in the ICCD photographs. The phenomena of stratification decreases 
in case of 100 mTorr due to the onset of shock wave ahead of the plasma front [78] which 
hinders the motion of the fast component. Furthermore, the dimension of the plume gets 
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Figure 4.5: The temporal profile of Cl transition (at 14 mm), CH transition (at 16 mm), 
Cm transition (at 14 mm) at 5 mTorr of air. The slow peak P 2 is dominant in all the 
three cases. 



Laser Produced Carbon Plasma in Low Pressure 


reduced, due to higher pressure from all sides and two separate fast and slow components 
are not resolved. However, in spectroscopy each species is observed individually so double 
peak structure is still observed in case of 100 mTorr whereas in case of ICCD photographs 
we observe the total emission from all the species, light emitted by fast and slow ions of 
different species merge together and the formation of stratification is lost. 

As the shock front expands in 100 mTorr air, more and more background gas is swppr 
up by the shock wave and the expansion velocity of the plasma plume decreases. This kind 
of non-steady shock wave is called a blast wave. Using Sedov’s approach, the motion of 


the blast wave front is governed by the relation Eqn.(1.24) 


where Co is a constant which depends on specific heat ratio of gas, E i is the laser energy 
and poj the density of the ambient gas. The above equation is derived for a strong explosion 
and assuming that the entire mass encompassed by the explosion wave is concentrated in 
a thin layer behind the front surface [50]. The R-t plot giving the displacement of the 
luminous front with time at 100 mTorr along with a fit to shock model S t = a 0 t p is shown 
in Fig.4.6. The value of p is found to be 0.28 comparable to the theoretical value of 0.4. 
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Figure 4.6: Experimentally observed plume front is fitted with the equation S t 




Chapter 5 


Laser Produced Carbon Plasma in a 
Nonuniform Magnetic Field 


5.1 Introduction 

The expansion of plasma produced by focussing an intense pulsed laser beam in a 
magnetic field is a rich and complicated phenomenum occuring in many laboratory envi- 
ronments. Am ong a number of external control parameters, magnetic field can be applied 
to control the local plasma thermofluid characteristics from both macroscopic anti mi- 
croscopic viewpoints [ 26 , 27 ] important for studying astrophysical processes and various 
technological applications 

This Chapter deals with the studies of LPP expansion in nonuniform magnetic field 
using three diagnostics viz. ICCD photography, time resolved and time integrated emission 

°The work presented in this Chapter is based on : 

1. "Laser-produced carbon plasma expanding in vacuum, low pressure ambient gas and nonuniform 
magnetic field". A. Neogi and R. K. Thareja, Phys. Plasmas 6, xxxx (1999). 

2. "Dynamics of laser produced carbon plasma expanding in a nonuniform magnetic field" A. N’eofi and 
R. K. Thareja, J. Appl. Phys. 85, xxxx, (1999). 

3. "Optical emission spectra of laser ablated carbon plasma in curved magnetic field" A. Neogi and R. 
K. Thareja, Communicated to J. Appl. Phys. 
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spectroscopy, and Langmuir probe. The Chapter starts with the basic theory of plasma 
behavior in magnetic field . The shape of the plasma expanding in vacuum, ambient gas 
and nonuniform magnetic field is observed through fast photography. The images are 
used to estimate the plasma parameters like electron density, electron temperature etc. 
[178]. The temporal profile of the different species of the carbon viz. Cl, CII. CIII CIV 
shows multiple peaks in nonuniform magnetic field [129], The time integrated .spectra of 
carbon plasma show the effect of curvature and gradient of magnetic field on the intensity 
of spectral lines and molecular C 2 emission [147]. Finally the use of Langmuir probe to 
monitor plasma behavior has been discussed. 


5.2 Laser Produced Plasma in Magnetic Field 

Energy, momentum and Ohm s law equation for the plasma pl um e expanding in magnetic 
field is given below [120] . 

5.2.1 Basic Equations of the Plasma in Magnetic Field 


(1) Energy equation for plasma in magnetic field is derived from first law of thermody- 
namics ; 


dQ_d£ dW 
dt dt dt 


( 5 . 1 ) 


where Q, U, and W have usual meaning defined in Chapter 3. and 

dQ 

dt 


Pl + W s 


( 5 . 2 ) 


where P L is the laser power and W 5 the heating due to surface current, and 


dU 

dt 


d_ 

dt 


■(N e + Ni)k B T 


3 , Ar dT 

= -(N e + N i )k B —, 


( 5 . 3 ) 


dW n dV {N' + NJksT A ,dr 

= P — = f x 47tr"— 

dt dt \nr z dt 


™ = l iN ' +Ni)T * 


or 


( 5 . 4 ) 



5. laser Produced Carbon Plasma in a Mammiform Magnetic 80 


and W, = J i 2 RdV = J ijJ 2 dV 


(5.5) 


where rj is the plasma resistivity given by Spitzer [193] and J is the surface current. 
Magnetic field is assumed to decay exponentially into the plasma with skin depth A, i.e. 
the depth at which field becomes 1/e of its original value (at the boundary) as it propagates 
a distance A; 

i.e B ~ Be for 0 < r*i < r 


(5.6) 


Since J = 


V xB 

47 r 


from Eqn.(5.6), it follows 




W s 


= J rjJ 2 dV = 


rjB 2 

87tA 


r 2 -rA + y(l-e~^) 


at ri = r 


(5.7) 

(5.8) 


Combining Eqns.(5.1), (5.2), (5.3), (5.4) and (5.8) we get the energy equation in final form 
as 

3 dT 3 dr nR 2 r . 1 

(5.9) 




o . A" , 2r . 

r~ — rX + —(1 — e T ) 


(2) The momentum equation is derived by equating the work done by the expanding plas- 
moid to the change in kinetic energy; 


dV A 

F S = 


(5.10) 


In the presence of magnetic field the above equation changes to 

P % -A(i^) + /( J xB)^± 


(5.11) 


where the second term in the R.H.S arises due to expansion in magnetic field. From 
Eqns.(5.6) and (5.7) it follows 

.dr r 9 „ > , m 2 /i g -£) 


J (J x B )dVj t = B 2 [r 2 - 2rA + 2A 2 (1 - e * 


dr 

dt 


( 5 . 12 ) 
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Combining Eqn.(5.4) and (5.11), (5.12), we get 

W + W| )l| - i -8±fe) . *| y _ 2rA + 2A , (1 _ e - {)] (5 13) 

(3) The Ohm’s law equation for the plasma motion in magnetic field is derived by equating 
the rate of change of total magnetic flux to the emf generated on the surface; 


dt 


(Flux) = - Jr) J.dl 


■qrB 

2~X 


(5.14) 


where dl is infinitesimal length along the plasma boundary 

Fiux = jf B(2nx)dx = jT 5e-^(27rx)dx = 2tBX [r - A(1 - e~$)] 


’ . Ohm’s law equation reduces to 


dt 


r 

47r A 


(5.15) 


Solving Eqn.(5.9), (5.13) and (5.15) for r, the radius of the plasma plume, one obtains an 
energy integral equation, 

yl'f-'-o] = §Wi + 


+2B 2 


I‘ F ^tA + jf w + (5 - 16) 


where 


m-r#i ” (W)I 

JO 


r 2 — rX + —(1 — e 2r / A ) 


3 87r(TV e + Ni)X(f) 
and F 2 (r , A) = r A - A 2 (l - e" r/A ) 


where ro and Tq are the initial radius and temperature of the plasmoid just after the laser 
pulse, N e /(Nj) is total number of electrons or ions, M the average mass of the plasmoid, 
B, the magnetic field and rt, is the radius of plasmoid where v = 0, this is usually referred 
to as "bounce" radius. The function fi(t) is proportional to the resistivity t) and F 2 (r, 
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A) depends on the skin depth A. The first and second terms on the ri^ht hand side of 
Eqn.(5.16) refer to thermal and kinetic energy respectively, third and fourth terms arise 
due to finite resistivity of the plasma. Radiation loses are not taken into account in arriving 
at Eqn. (5.16). Solving for r h , we have 

r t 3 ~ i^f + | (N, + Ni)k B n+HH |F,(r.A)(it 

+-p(K+N i )k B &l#L$dt (5.17) 

It follows from Eqn.(5.17), the plot of r h with time has an oscillatory nature. The radius 
of the plasma plume shows periodic bounce, i.e positions of v = 0. If the resistivity is 
zero the bounce radius remains constant and lines of force cannot diffuse into the plasma 
whereas if the resistivity is finite the magnetic field line slowly diffuses into the plasma and 
the bounce radius increases with time and oscillations stops. 

5.2.2 Diamagnetic Behavior in Magnetic Field 

The motion of plasma in magnetic field is analogous to a number of loops moving through 
a magnetic field. For a loop travelling through a magnetic field, the induced emf £ is .given 

by 

£ = (5.18) 

at 

where <b is the flux, I, the current and R is the resistance. Plasma being conducting fluid, 
we can assume, R = 0. Eqn. (5.18) gives 



or 


(p = Constant 


(5.19) 
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It implies that magnetic- flux through any surface in the fluid tends to remain constant 
as that surface moves with the fluid. Lines offeree are "frozen in" with conducting fluid. 
Motion of field lines with respect to the plasma is possible only if resistivity is not zero. 
The magnetic' flux across the plasma surface can change only if the magnetic field diffuses 


into the plasma winch is determined by the resistivity of the plasma. The magnetic field 
behaves diamagnetically i.e. the field within the plasma is less than that in the surrounding 
space. This results in the retardation of the motion of the plasma across the field lines. 


5.2.3 () and Magnetic Reynold’s Number R m 

As the plasma moves in a magnetic field B, the charged particles experience a force per- 
pendicular to the motion of the plasma and applied magnetic field due to which a surface 
current is generated perpendicular to the moving plasma. Mathematically, 

F ■= nqv x B (5.20) 


where n is the number of particles 
or 

F=JxB (5.21) 


where J = nq v, is the surface current. 

From Maxwell equation 

V x B 
— r 

4ir 

JxB = B X V x B 

4-7T 


But 


jd2 

B x V x B = V(— •) - (B.V)B 

jL 


On substituting Eqn.(5.24) in Eqn.(5.23) we get 

J X B = -V (fi) + i(B.V)B 
8-7T 47T 


{5.22) 

(5.23) 


(5.24) 


(5.25) 
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w ^ erP S 7 r * s ma K not:ic energy density which, plays the role of magnetic pressure. The 
term £(B.V)B .irises due to curvature of the field. If the field is unidirectional i.e, B . V 
= 0, Eqn.(5.26) reduces to 

1>B = ~ V <S> (5.26) 

3 is defined as the ratio of the particle pressure to the magnetic pressure, viz. 


(3 = 


ParticLcPressure 


nkaT 


(5.27) 


MagneticP res sure B 2 /8n B 2 /8 tt 

In our experiments we have. n~ 10 17 / cm 3 , T~ 7eV and B = 3 kG at 2 mm from the 
target (section 4.2) and hence fi comes out to be 3. The plasma pressure being higher than 
the magnetic pressure, the plasma "punches" its way across the magnetic field. As the 
plasma expands the (3 value rapidly decreases, approaching unity. Magnetic field simp ly 
diffuses into the plasma and the plasma becomes electrically polarized, the electric field 
E is established which is perpenduicular to both B and the plasma velocity. The plasma 
moves across the B field by E x B drift. The plasma expansion ceases when the displaced 
magnetic energy approximately equals the kinetic energy of the plasma expansion. 

Magnetic Reynolds number R m determines the degree of magnetic field distortion of 
lines. For Rm > 1, the background field lines ahead of the flowing plasma will be distorted 
and compressed whereas for Rm < 1, the external field lines will be relatively unaffected 
and simply diffuse through the expanding plasma. 

To get Rm we start with the equation for the evolution of the magnetic field 

= V x (v x B) + n m V 2 B 

at 


(5.28) 


where is called the magnetic viscosity defined by 

1 

Um (4ir ctq) 


(5.29) 


where <r 0 is the longitudinal conductivity i.e. conductivity parallel to the B field [6]. The 
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first term on the right hand side of Eqn.(5.28) is railed the flow term while the second 
term is called the diffusive term. The ratio of the flow term to the diffusive term gives the 
magnetic Reynolds number [6], viz. 


where v is the velocity of the fluid plasma, L the magnetic field gradient scale length and 
n m the magnetic viscosity. 

Using our experimental value we get 5 i.e the flow term dominates over the 
diffusive term which inturn results in compression of magnetic field lines. 


5.3 Experimental Techniques 

The experimental set up used in the present study is described in Chapter 2. The field is 
maximum (3.5 kG) at the middle of the poles and decreases on either sides (along X axis) 
whereas along the Z direction it is minimum at the middle and ma ximum (5.1 kG) at the 
pole surface. The distance between the poles is 3.2 cm. The target is kept in between the 
pole pieces at one extreme end. The energy of the laser pulse used is 800 mJ /pulse having 
wavelength of 1.06 n m. 


5.4 Results and Discussions 

5.4.1 ICCD Photography 

The plasma exp ansi on in vacuum (< 1 mTorr), 100 mTorr air and non-uniform magnetic 
field in vacuum was recorded by ICCD camera at different delay times. Figure 5.1 shows 
ICCD images at different delay times. The exposure time for all these photographs is kept 
constant (20 ns) so that relative intensity can be compared. In order to view the real size 
of the images of varying intensity, each plume is normalized with respect to its maximum 
intensity. 

The displacement of the luminous plasma front with respect to the target position 
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Figure 5.1: ICCD pictures of the plasma at different delay times with 20 nsec exposure. 
The images are normalized with respect to the ma ximum intensity of each image, fa) 
Vacuum, (b) 100 mTorr air and (c) Variable magnetic field in vacuum. 












100 1000 


Time ( ns ) 

Figure 5.2: The variation of the plume front with time in vacuum, air at 100 mTorr and 
magnetic field. The dotted lines connecting the points are only an aid to the eye. 

is plotted as a function of time in Fig. 5. 2. The volume of the expanding plasma plume 
for three different conditions is shown in Fig. 5. 3. The error in these quantities mainly 
originate from the error in determining the exact location of the plasma boundary. 

The shape of the plume in vacuum is spherical throughout. The plume expansion is 
highly forward directed in this case. The plume appears very much diffused after 500 us 
with 20 ns exposure time as shown in Fig.5.1(a). At 100 mTorr air, the shape of the 
plasma plume remains semispherical but the expansion of the plume gets highly reduced 
due to larger thrust from all sides as shown in Fig.5.1(b). The temporal behavior of plasma 
plume in the presence of nonuniform magnetic field in vacuum is shown in Fig.5.1(c). The 
size of the plasma as it propagates shows alternate contraction and expansion due to which 
the temporal behavior of the plasma front shows oscillations as shown in Fig.5.2. The first 
compression is seen at 150 nsec and the second at 500 nsec, Fig.5.1(c). To understand 
alternate contraction and expansion consider Eqn.(5.16). If we neglect resistive effects and 
assume that at early time y Vq << f (N e + Nj) ks To and r(, >> ro, we get 
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Figure 5.3: Temporal variation of the volume of the plasma plume in vacuum, air at 100 
mTorr and magnetic field. The lines connecting the points are only an aid to the eye. 
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The period of oscillation being given by 

r b 

T ~ — . 

V 


(5.31) 


(5.32) 


where v is the velocity of the plasma front. In order to calculate Vb and r [Eqn.(5.34) 
and Eqn.(5.35)], it is assumed that T 0 during initial expansion is not influenced by the 


magnetic field [118, 203] and hence 


T 0 < 


M \v 2 

5 ( 1 + 1 ) 


(5.33) 


where M, = 12 x 1.67 x 10 -27 kg, v = 1.1 x 10 5 m/sec is asymptotic plasma expansion 
velocity measured from the ICCD photographs, taking z = 2, T 0 is estimated to be 101 
eV from Eq(36). Taking N = 5.48 x 10 16 (calculated from the ICCD photographs) and B 
= 3 kG we get v h « 1.64 cm and r = 0.04 p sec. From Fig.5.2 which shows the variation 
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of position of the plume front at various time, we find r b to be 0.43 cm and r to be 150 
ns, different from those calculated using Eqn.(5.31) and Eqn.(5.32). The radius of the 
second bounce is smaller ( Fig.5.1. Fig. 5. 2 ] and the time period larger. The decrease in 
bounce • radius with time in our ease may be due to the energy loss from the plasmoid 
due to Bremsstrahlung radiation which is quite high in the early phase of expansion and it 
dominates over the resistive terms. This effect is not accounted for in Eqn.(5.16). Secondly, 
magnetohydrodynamic instabilities such as Rayleigh-Taylor instability which damps the 
plasma motion are also not considered in the equations of the plasma dynamics. 

It is observed in Fig.5.1(d) that at delay times greater than 800 ns, the plasma moves 
towards one of the poles. The ICCD photographs of the plasma taken at a delay of 1000 
ns when the target is in the middle and shifted towards the upper pole of the magnet is 
given in Fig.5.4. It is clear from the photographs that at a later time, the plume breaks 
into two lobes which move towards the poles. Figure 5.1(d) corresponds to the target 
position 3 nun nearer to the lower pole. When the target position is in between the poles, 
the lobes are symmetric, otherwise asymmetric. The above fact can be explained using 
magnetohydrodynamic model. 

The generalised form of Ohm’s law in a magnetoionic fluid is given by [6] 

E + v x B = J /(To + (J x B )/n e e (5.34) 

where E, B, J, . < 7 () , n,,, and v have usual meaning defined earlier. In an one dimensional 
flow, in the absence of conduction currents E y = v^EL. As the plasma plume expands in 
transverse magnetic field, the ions and electrons get separated due to Lorentz force till the 
magnetic force acting on the charged particles equals the electric force. The plasma moves 
across crossed electric and magnetic field with a drift velocity given by 

v = E x .— (5.35) 

B 2 

On account of three dimensional expansion of the plume, the induced electric field trans- 
verse to the flow will vary over the plume cross section because v z decreases at larger 
angles. The field gradient in the plume boundary allows the closed current paths to form. 
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Figure 5.4: The expansion of the plasma in three different positions of the target at 1000 
ns delay with respect to the ablating laser pulse, (a) Target is at the middle of the poles; 
(b) the target is few mms towards the upper pole of the magnet; (c) target is more towards 
the lower pole. 
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Figure 5-5: A cross sectional area of the plume showing the direction of current within 
and at the boundary. The J x B force at the boundary accelerates the plume boundary 
giving rise to lobes which move towards the poles. 

The current which flows inside the plume interacts with the magnetic field giving rise to 
J x B term which decelerates the flow. Outer regions of the plume would likewise be 
expected to accelerate as return current flows in the opposite direction leading to -J x 
B term. This J x B term increases with increase in B with larger Z (towards the poles 
in Fig. 5. 5) so the plume is observed to break into two lobes which move towards higher 
magnetic field side. The flow of induced current and formation of lobes is schematically 
shown in Fig. 5. 5. 

The experiment is also done in magnetic field at 5 mTorr and 100 mTorr air pressure. 
In both the cases an increase in intensity is observed in addition to the features observed 
without the magnetic field. 

The velocity profile of the plasma front in vacuum, 100 mTorr air and magnetic field 
is shown in Fig. 5. 6. The plume expansion is highly forward directed in case of vacuum 
with leading edge moving with velocity of 1.1 x 10' cm/sec at 29 ns after the ablating 
laser pulse. The expansion velocity is highest in vacuum followed by air at 100 mTorr 
and magnetic field. Oscillations are observed in case of magnetic field due to alternate 




100 1000 
Time ( ns ) 

Figure 5.6: The velocity of the plume front with time. 


contraction and expansion as discussed earlier. Velocity in magnetic field is minimum ac 
150 nsec and 500 nsec when the radius of the plasmoid is lowest. 

Variation of mean and maximum intensity of the l umi nous plume is shown as a func- 
tion of time in Fig. 5. 7 (a) and (b). The temporal behavior of I memi and 1^ are similar. 
The intensity in case of 100 mTorr air is highest followed by magnetic field and vacuum. 
In case of magnetic field a sudden rise in intensity is seen at the bounce points (150 ns, 500 
ns) because of contraction of the plasmoid. At 900 nsec an increase in emission intensity 
is observed but no contraction is observed. The small fluctuations in intensity for the case 
of magnetic field in the interval 300-1200 nsec is shown in the inset of Fig.5.7 (b). 

The ICCD images were used to estimate the electron density and temperature of the 
plasma. Referring to hydrodynamic Eqn.(l.ll, 1.16) before and after the te rmina tion of 
the laser pulse, we find that after the termination of laser pulse, expansion in the direction 
of propagation that is in X direction, is well beyond the expansion in Y and Z direction. 
Since the laser pulse is Gaussian, the expansion in Y and Z direction can be taken as 
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Figure 5.7: The t em poral behavior of maximum and mean intensity of ICCD photographs. 
The inset in fig. (b) shows expanded version of the variation of intensity in case of magnetic 
field. 
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Figure 5.8: The temporal variation of electron density calculated from the ICCD pho- 
tographs. 

symmetric. From mass continuity equation it follows, 

p{to)V(t 0 ) = p(t)V(t) (5.36) 

Here p is the ion density , V, the volume of the plasmoid, to (< t)refers to time for the 
photograph taken at the earliest time and t is any later time. Electron density at an 
early time is calculated theoretically assuming adiabatic expansion [107] and volume at 
early time is calculated from the photographs. Thus knowing the volume of the plume at 
any time (from the photographs), electron density can be estimated. Figure 5.8. shows 
the temporal behavior of calculated electron density. The temporal behavior shows thac 
electron density is highest in the presence of magnetic field followed by air and vacuum. 
Further, it shows oscillatory behavior in magnetic field due to contraction and expansion 
of the plasmoid. 

The temperature of the plasma was also calculated from ICCD photographs. Since 
the temperature is high it is assumed that the plasma acts like an assembly of black body 
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radiators and the emission is mostly in the visible region (due to high temperature, Wiens 
displacement law). According to Stefans’law, the power radiated (P r ) by a black body is 

a T 4 S = P r ex total intensity observed in the photographs (h) 

Thus at t and t 0 one can write 

crT 4 (t)S(t ) = KiJ T (t), t > t 0 . (5.37) 


<rr 4 (t„)S(t 0 ) = KMk), t = f„ (5.38) 

Here, S is the total surface area and K b the proportionality constant. The temperature at 
early time is calculated theoretically assuming adiabatic spherical expansion at early time 
[107] and area is calculated from the photographs. Comparing Eqn.(5.37) and Eqn.(5.38), 
it follows 


T 4 (t)S(t) T\t 0 )S(t 0 ) 

~m~ = TT ( ' U9) 

Calculating the surface area at various times from the photographs, the temperature was 
calculated and is shown in Fig.5.9. It is found that the temperature in 100 mTorr air 
is highest followed by magnetic field and vacuum. A slight increase in temperature in 
magnetic field is observed at the second bounce point (500 ns). Increase in temperature in 
gas is caused due to enhanced interparticle collisions due to compression from the outside. 
In case of magnetic field, surface currents are set up due to finite resistivity at the plasma 
field boundary producing increased temperature on account of ohmic heating. The extra 
energy due to o hmi c heating comes from the plume kinetic energy. 

Assuming adiabatic expansion of the plasma plume after the ter m i n ation of laser pulse, 
the theoretical calculations show that electron density and electron temperature in vac uum 
vary with time as meet -3 and T e at -2 [107]. We have fitted the equation, y — w ^ r ^- 

the calculated values of electron density and electron temperature for the vacuum case as 
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Figure 5.9: The temporal behavior of electron temperature calculated from the ICCD 
photographs. 

shown in Fig.5.10. The value of p comes out to be -2.85 and -2.73 close to the theoretical 
value of -3 and -2 for the electron density and electron temperature respectively. 

5.4.2 Time Resolved Emission Spectroscopy 

The emission spectra of carbon plasma were recorded at different distances away and 
parallel to the target surface by moving the monochromator in a plane perpendicular to 
the target surface. The temporal profiles of Cl transition 2p 3 3 D° - 8f F (5/2) at 399.7 nm. 
CII transition 3d 2 D - 4f 2 F° at 426.7 nm, CIII transition 3s 3 S - 3p 3 P° at 465.0 nm and 
CFV transition 3s 2 S - 3p 2 P° at 580.1 nm were recorded at different distances away and 
parallel to the target surface in vacuum (P = 0.1 mTorr), with and without magnetic field. 
The emission lines were identified using the information available in the literature [186;. 

In the absence of magnetic field the temporal profile of all the species showed a single 
peak structure as shown in the inset of Fig. 5. 11 - 5.14. In the present case, the plasma 
moved in a non-uniform magnetic field as shown in the inset of Fig.2.1. (Chapter 2). Along 
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Figure 5.10: The n ( . and T e of the plasmoid expanding in vacuum, calculated from ICCD 
photographs has been fitted with equation of the form, y = b 0 t p . 

Z direction, the field was maximum at the poles and mi nimum in the middle. The single 
plume was observed to break into two symmetric lobes (say lobe 1 or Ll and lobe 2 or 
L2), one towards the north pole and other towards the south pole as shown in Fig.5. 15(a). 
However, it was very difficult to analyze the different components (fast/slow), arising from 
each of these lobes spectroscopically because the two lobes (Ll. L2) being symmetric, fail 
on the same line on the detector. To see the effect, the target was kept off centered such 
that the laser was focused nearer (~ 3 mm) to one of the poles than the other which gave 
rise to asymmetric lobes, Ll and L2, Ll being larger than L2. Figure 5.15 (a) and (b) 
(at the left) show the observed plume at 1000 ns with target kept at the centre and off 
centred with respect to the poles respectively. The Fig.5. 15(c) and (d) give a pictorial 
representation. 

The temporal profiles of CII and CIV ions which are at the front of the plasma plume 
show interesting features at a distance of 6 mm onwards from the target surface. Similar 
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Figure 5.11: Time resolved emission spectra of CII species in vacuum with magnetic field 
at (a) 4 mm, (b) 6 mm, (c) 8 mm and (d) 10 mm respectively along with the best fitted 
curve (solid line) and the fitting Gaussian functions (dotted lines) .In the inset the dotted 
line represents the emission spectra without magnetic field and solid fine represents FFT 
smoothed spectra of the magnetic field case. 
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Figure 5.12: Time resolved emission spectra of CIV species in vacuum with magnetic field 
at (a) 4 mm, (b) 6 mm, (c) 8 tutu and (d) 10 mm respectively along with the best fitted 
curve (solid line) and the fitting Gaussian functions (dotted lines).In the inset the dotted 
line represents the emission spectra without magnetic field and solid line represents FFT 
smoothed spectra of the magnetic field case. 
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Figure 5.13: Time resolved emission spectra of CIII species in vacuum with magnetic field 
at (a) 14 mm, (b) 16 mm, (c) 18 mm and (d) 20 mm respectively along with the best fitted 
curve (solid line) and the fitting Gaussian functions (dotted lines). In the inset the dotted 
line represents the emission spectra without magnetic field and solid line represents FFT 
smoothed spectra of the magnetic field case. 
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Figure 5.14: Time resolved emission spectra of Cl species in vacuum with magnetic field 
at (a) 14 mm, (b) 16 mm, (c) 18 mm and (d) 20 mm respectively along with the best fitted 
curve (solid line) and the fitting Gaussian functions (dotted lines). In the inset the dotted 
line represents the emission spectra without magnetic field and solid line represents FFT 
smoothed spectra of the magnetic field case. 
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Figure 5.15: (a) ICCD images of the plume breaking into two symmetric lobes (at 1Q0U nsj 
when focal spot is exactly between the two poles, (b) Plume breaking into two asymmetric 
lobes (at 1000 ns) when focal spot is nearer to one of the poles, (c) Pictorial representation 
of (a), (d) Pictorial representation of (b). LI , L2, C refers to lobe 1, lobe 2 and central 
lobe. F1/(F2), 11/(12), S1/(S2), refers to fast, intermediate and slow components of lobe 
1 and lobe 2 respectively. 
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features are also observed in CIII and Cl at 16 mm and onwards. The time resolved spectra 
of CII, CIV, CIII and Cl ions along with the best fitted curve at different distances in the 
direction of propagation of the plasma plume are shown in Fig.5.11 - 5.14. respectively. 
It has been found that Gaussian function fits best with the different peaks observed in 
the spectra. The best fitted curve (solid line) (minimum standard fit error of around 
0.01) is the resultant of individual Gaussian curves (dotted line) shown in the figure. The 
significance of fitting the main spectra with these individual curves lies in the fact that the 
peak position of these individual curves help locate the less prominent peaks in the original 
spectra in addition to the prominent ones. It also helps understand their evolution with 
time and distance. Figure 5.11 shows time resolved spectra of CII at 4 mm and onwards up 
to 10 mm along with the best fitted curve. In the inset the solid line shows the smoothed 
curve by feist fourier transform (FFT) filtering, and dotted line shows the temporal profile 
in vacuum without magnetic field plotted on the same scale as that of magnetic field for 
comparison. It shows that intensity of the spectra increases in presence of magnetic field. 
After 10 mm , the intensity of all the peaks decreases till 14 mm beyond which signal is 
too weak to be recorded. The spectrum of CIV from 4 mm onwards is shown in Fig.5.12. 
The spectrum of CIV like that of CII could be recorded till 14 mm. The spectrum of CHI 
at different distances is shown in Fig.5.13. Till 14 mm. only one single peak is observed 
which can be fitted by 3 Gaussian functions like that of CII at 4 mm. At 16 mm that is 
just after the disappearance of CII and CIV, CIII shows small peaks in addition to the 
main peak. The case is very much similar to that of CII and CIV. The spectrum of Cl 
shown in Fig. 5. 14 resembles CIII very closely. 

The velocity of the carbon species is calculated from the delay of the peak of the 
emission spectra at different distances with respect to the trigger pulse. The velocity of 
the carbon species with and without magnetic field at a distance of 4 mm from the target 
is shown in table 5.1. In magnetic field, velocity of CII species is highest followed by CIV, 
CHI and Cl species. 

In addition to the multiple peak structure, another interesting feature observed is the 
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Table 5.1: Velocity of carbon species with and without magnetic field in vacuum. 


Species 

Velocity (10 6 cm/sec) 


B absent 

B present 

Cl 

1.58 

1.54 

CII 

1.24 

2.50 

CIII 

1.76 

1.68 

CIV 

1.90 

1.69 


oscillations in the temporal profile of CII and CIV at distance of 6 mm to 14 mm away from 
the target and at 16 mm to 20 mm for CIII and Cl respectively. Since these oscillations 
are observed for the species at the outermost boundary they can be attributed due to edge 
oscillations in magnetic field. 

The formation of multiple peaks can be explained using magnetohydrodynamics model 
as discussed earlier in section 5.41 [Eqn. (5.34)] On account of three dimensional expansion 
of the plume the induced electric field transverse to the flow will vary over the plume 
cross section because v x decreases at larger angles. Field gradient will exist in the plume 
boundary permitting closed current paths to form. The current which flows inside the 
plume interacts with the magnetic field giving rise to J x B term which decelerates the 
flow. This gives rise to slow component (S) in the plasma plume. In the outermost 
boundary, the electrons flow in the opposite direction giving rise to -J x B term which 
accelerates the outermost boundary [30] giving rise to fast component (F) in plasma plume. 
The region lying in between fast and slow component forms the intermediate component 
(I). It is schematically shown in Fig.5.16. 

Multiple peak formation in the emission spectra is due to two lobes, LI and L2, which 
move towards the poles and the central lobe. Each lobe has three components viz. the fast, 
the intermediate and the slow one. In Fig.5.11 at 4 mm from the target, other than the 
first small peak F0, there is only one distinct peak which can be best fitted ( least standard 
fit error of around 0.01) by three Gaussian peaks shown as F (fast), I (intermediate) and S 
(slow) . The s mall peak F0 is the one also observed without the magnetic field. On going 
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Figure 5.16: Schematic diagram of the direction of current density J giving rise to fast, 
intermediate and slow component. 

from 4 mm to 10 mm, it is observed that the fast peak (F) breaks into three distinct peaks 
(FI, F2, F3), the intermediate (I) and slow peak (S) break into two distinct peaks (II, 
12) and (SI, S2) respectively. FI, II, SI are the fast, intermediate and slow components 
coming from LI whereas F2, 12, S2 are due to L2. F3 is the fast peak due to central 
lobe. The intermediate and slow peaks of central lobe are not resolved. The situation is 
schematically depicted in Fig.5.16. In general we observe 

(1) Peaks due to LI are of higher intensity than the corresponding peaks from L2 i.e. 
intensity of F1>F2, Il>12 and Sl>S2. This is because intensity of LI is higher than that 
of L2. (Fig.5. 15(b)). 

(2) The intensity of the fast peak F3 due to central lobe is of lower intensity than F1/(F2) 
and the intensity decreases faster compared to F1/(F2). This is because central lobe 
disintegrates into LI and L2 with increasing distance. 

Intensity and velocity are calculated for all the peaks of each of these species and their 
nature is studied. The velocity of the fast peaks for CII (Fig.5. 11), CHI (Fig.5. 13) and 
Cl (Fig.5. 14) increases with distance whereas for other peaks, it decreases with distance 
For CIV species velocity increases for all the components but it is highest for the fast 
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Figure 5.17: Variation of electron temperature with distance with and without magnetic- 
field. 

components followed by intermediate and slow components. 

The effect of magnetic field on CIH and Cl species is conspiquous only after the 
disappearance of CH and CIV species after 14 mm when CIH and Cl species are at the 
outermost boundary. In both CHI ( Fig.5.13) and Cl (Fig.5.14) till 14 mm, only one peak 
is observed in the signal which can be best fitted by three Gaussian peaks. However. 16 
mm onwards, multiple peaks are observed in CIII and Cl, similar to that of CII at earlier 
distances. Beyond 20 mm , the spectrum of CIH and Cl could not be recorded. 

The increase in intensity of the temporal profile of all these species ( Cl, CH, CIH 
and CIV) in magnetic field is due to higher temperature and density of the plasma in 
magnetic field. The temperature profile of plasma with and without magnetic field is 
shown in Fig.5.17. The temperature of the plasma was experimentally determined at 
various distances using relative line intensities of the species having same ionization state 
assuming the plasma to be in local thermodynamic equilibrium. The temperature of the 
plasma is found to be higher in magnetic field case. The increase in electron temperature 
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may be due to joule heating and electromagnetic compression of the plasma. The density 
can be estimated from model of mass ablation rates [133], that is dm/dt a l£ 6 . Since 
the velocity of the bulk of the plasma decreases in magnetic field, the density increases in 
magnetic field compared to no magnetic field case. The electron density at 4 mm from the 
target is estimated to be 1.18xl0 17 / cm 3 without magnetic field compared to 1.24 xlO 17 / 
cm 3 in magnetic field. Increase in density in magnetic field has been reported earlier [177] 
The temporal variation of density and temperature as shown in Fig. (5.8, 5.9) shows rise 
in both density and temperature in magnetic field. 

Oscillations are observed in the temporal profile of CII and CIV from 6 mm and 
onwards (Fig.5.11 - 5.12); in CHI and Cl from 16 mm and onwards (Fig.5.13 - 5.14). These 
oscillations are not observed without magnetic field. Since the oscillations are observed for 
species when they are at the outermost front, they are attributed due to edge instability 
in magnetic field. These edge oscillations are similar to that observed by Kent eta i [136] 
in a transverse magnetic field of 4.4 kG using langmuir probe. 

5.4.3 Time Integrated Emission 

T ime integrated spectra of carbon plasma is recorded at various distances from the target 
surface kept at positions marked as P, Q, R respectively between the poles, as shown in 
Fig.5.18. The position of these points and the respective magnetic fields at these points is 
also shown in the figure. Figure 5.19 shows the time integrated spectra of carbon plasma 
recorded at 6 mm (P'/ Q'/ R') from the target position (P/ Q/ R) respectively in the 
spectral range 350 nm to 600 nm, thus the spectra recorded at position P' correspond to 
target position at P (Fig.5.18) and similarly for other cases. The spectra without magnetic 
field are also shown in Fig.5. 19(d) for comparison. The wavelength and the corresponding 
excited ionic states observed are shown in the figure. The emission lines are identified using 
the information available in the literature [184]. The characteristics of emission spectra 
vary considerably for the three positions in the magnetic field and that of without magnetic 
field. The striking differences observed in the four cases are enumerated below. 
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Figure 5.19: Time integrated emission spectra at three different positions in magnetic field 
(a) target kept at P and spectrum taken at P'; (b) target kept at Q and spectrum taken 
at Q'; (c) target kept at R and spectrum taken at R' and (d) B = 0 case. 
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The intensity of the emission lines increases in magnetic field for the position P and Q 
of the taxget compared to that of no magnetic field case. It is highest in case of position P 
followed by position Q of the target in magnetic field. However, the intensity of the lines 
in position R of the target is less than that of no magnetic field case. 

The excited lines are both from neutrals and ionic species but the concentration of 
ionic species is more than that of the neutrals. The intensity variation of a strong line of 
CII (426.7 run, 3d 2 D - 4f 2 F°) and CIII (407.0 nm, 4f 3 F° - 5g 3 G) is plotted as a function 
of distance for all the three position of magnetic field in Fig.5.20. It follows from the figure 
that the intensity of the fines increases initially, attains a maximum value and then starts 
decreasing. The fall in intensity with distance for the spectra recorded at position P is 
slower than that for position R.The emission was observed as far as 30 mm from the target 
for position P compared to 10 mm for position R. 

It is interesting to note that the intensity of C 2 Swan band emission, d 3 7r ff - a 3 7T, x , with 
A v = -1, 0, 1 is much larger at the position P of the target compared to other positions 
as shown in Fig.5.19 and Fig.5.21. Figure 5.21 shows the band emission along with the 
transitions on an enlarged scale for A v = -1. 0, 1. Variation of C 2 band intensity (->16.5 
nm d 3 ? Tg - a 3 7 r u ) with distance for various positions and that without magnetic field is 
shown in Fig.5.22. The C 2 band intensity is highest in position P followed by position 
Q. The vibrational temperature of C 2 molecule was calculated from intensity ratio of two 
transitions [205] and was found to be 1.96 eV at a distance of 6 mm from position P. 

The temperature of the plasma is calculated from the ratio of relative intensity of 
the lines of the ionic species assuming plasma to be in local thermodynamic equilibrium 
(LTE) . For LTE to hold, the electron density must satisfy [45] 

n e > 10 n T; /2 (E k - E^cm^ (5.40) 

where T e is in k and (E^ - E,) in eV. 

Using A E = 2.84 eV, corresponding to the largest energy gap of transitions used for 
temperature measurement and highest T e (11.4 eV) the lower limit for n e from above 



5. Laser Produced Carbon Plasma in a Nommiform Magnetic. 


Ill 



Distance (mm) 



55 

c 

CD 

«*■* 

c 



Figure 5.20: Variation of relative intensity of the lines of (a) CII (426.7 nm) and (b) CIII 
(407.0 nm) transition with distance when target is kept at different positions in magnetic 
field and without magnetic field. 
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Figure 5.21: C 2 Swan, band emission at 6 mm from the target when it is kept at position 
P of magnetic field. Transitions corresponding to (a) A v = -1; (b) A v = 0 and ( 0 ) A v 
= 1 are shown in an enlarged scale. 
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Figure 5.22: Variation of relative intensity of C 2 band (516.5 am) with distance when 
target is kept at different positions in magnetic field and without magnetic field. 

equation comes out to be 8.33xl0 lj / cm 3 . Our observed values of n P calculated from the 
Stark broadened profile of CII transition is always greater than this limit which justifies 
the use of LTE approximation for our analysis. Figure 5.23 shows the spatial variation of 
temperature in presence and absence of magnetic field. The temperature in case of position 
P is found to be highest followed by position Q and R. It is estimated to be 11.4 eV, 5.10 eV 
and 4.34 eV at a distance of 6 mm from the target (P ; , Q', R') when the target position is P 
, Q and R of magnetic field respectively. The temperature in the absence of magnetic field 
at 6 mm is estimated to be 4.5 eV. Temperature is also calculated at different laser fiuence 
in magnetic field (in position Q) and without magnetic field at a fixed distance from the 
target (4 mm) as shown in the Fig.5.24. It is found that temperature is higher in magnetic- 
field and increases with increase in fiuence. The estimated density using Stark broadened 
profile without magnetic field at 4 mm from the target is 1.18 x 10 17 /cm 3 compared to 
1.24 x 10 17 /cm 3 at the same distance in the position Q of the magnetic field in agreement 
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Figure 5.24: Variation of electron temperature with laser fluence with and without mag- 
netic field at a fixed distance (4 mm) from the target surface. 

with earlier report [177], It is due to the compression of the plasma by JxB force and 
enhanced ionization due to rise in temperature in magnetic field. 

The visible observation of the colour of the luminous plume revealed that it changes 
from blue to green in the presence of magnetic field. The C-_> emission at 516.5 am may be 
responsible for this green colour in magnetic field. As is obvious from Fig.5.19 that in the 
absence of magnetic field, most of the strong emission lines fail in the wavelength range 
360 run to 465 nm which is violet and blue region in the electromagnetic spectrum. 

It is worthwhile to mention that in case of position Q of the target, the plasma traverses 
a path of highest magnetic field (QQ') followed by position R and P (Fig.5.18) whereas 
intensity is highest in position P followed by position Q and R. This observation signifies 
that curvature and gradient of magnetic field is responsible for the above observations. 
Field lines axe convex to the moving plasma with positive field gradient in position P 
compared to position R where field lines are concave to the moving plasma with negative 
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Hpld gradient. The earlier experiment in magnetic field [144] corresponds to position Q of 
target iu our experiment. 

In order to explain the characteristic emission in magnetic field we believe the following 
three processes may be responsible viz, 

(1) The compression of magnetic field due to flowing plasma which is different for three 
positions P, Q and R. 

(2) Flow of current due to electro-magnetic interaction. 

(3) Drift of charged species due to curvature and gradient of magnetic field. 

The expanding plasma has finite conductivity that can support currents, the plasma 
hydrodynamic motion leads to compression and distortion of the ambient magnetic field. 
For the plasma having high d, that, is when the particle pressure is higher than magnetic- 
pressure compression of field lines is more 1 . A better way of quantifying this fact Is by 
magnetic Reynolds number, Ify which is the ratio of magnetic field convection to magnetic 
field diffusion given in Eqn.(5.30). 



n m 


where v is the fluid velocity and L. the magnetic field gradient scale length and n,„, the 
magnetic viscosity. In our case at early time v ~ 10 a m/sec, L ~ 0.6 cm and n^ = 1.35 x 
10 2 m 2 /sec, we get R™ ~ 5 which is >1, that is the flow term dominates over the diffusion 
term which results in compression of the magnetic field lines. Figure 5.25. schematically 
shows the distortion of magnetic field lines for three positions. The field lines are highly 
compressed at position P where radius of curvature is opposite to the moving plasma 
whereas they are rarefied or less compressed at position R where the field lines are along 
the moving plasma. Thus due to distortion of field lines the effective magnetic field becomes 
highest at P followed by position Q and R. 

The flow of current due to electromagnetic interaction is governed by generalised 
form of Ohm's law in a magnetoionic fluid given by Eqn.(5.16) As the plume expands 
three dim ensio nall y, electric field induced transverse to the flow vary over the plume cross 
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Figure 5.25: Sch em atic diagram of the distortion and compression of magnetic field lines 
when target is kept at different positions in magnetic field. 


section, field gradients exist in the plume boundary and current starts to flow. 

In a curved vacuum field, the drift experienced by the charged species due to curvature 

and gradient of magnetic field is given by [202]. 

m Rcr x B 

V d -V C + Vg = 


q RLB* 


,1 - 
v + 


( 5 . 41 ) 


where v c and v g are the drift due to curvature and gradient of the magnetic field. v,,and 
v denotes velocity parallel and perpendicular to magnetic fines of force, the radius of 
curvature and B is the magnetic field. 

As discussed above the behavior of the charged species in magnetic field is mainly 
governed by two factors viz. electromagnetic interaction [Eqn.(5.16)j and curvature gradi- 
ent drift [Eqn.(5.41)] along with magnetic field distortion [Eqn.(5.28)j. Since mass of the 
electron is much less than that of ions, movement of ions due to electromagnetic force (acc 
= qvB /m) can be neglected compared to that of electrons. Similarly, drift for the electrons 
can be neglected compared to that of ions due to curvature of magnetic field (y d x m). 
In position P, the electrons and ions move in the same direction compared to position R 
where they move in opposite direction because of the reversal in the direction of radius of 
curvature. The situation is schematically shown in Fig.5.26. The movement of electrons 
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Figure 5.26: Motion of electrons and ions due to electromagnetic interaction and curvature 
gradient drift (a) positive field gradient with radius of curvature opposite to the direction 
of propagating plasma; (b) with the gradient less than case (a); (c) negative field gradient 
with the radius of curvature in the same direction of propagating plasma and (d) no 
magnetic field. 
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and ions in the same direction may favor recombination of electrons with ions to form C > 
molecule in position P compared to position R. This may be the reason for not getting C_. 
emission at position R. Regarding the. motion of electrons and ions, position Q is similar 
to position P but the velocity of the electrons is less in position Q compared to position P 
because of lower effective magnetic field in position Q. This results in generation of less C, 
in position Q compared to position P. The velocity of the electrons is highest in position P 
followed by position Q and R because of difference in effective magnetic field which causes 
ohmic heating highest in position P resulting in highest intensity of spectral lines followed 
by position Q and R. The reason for lower intensity in position R compared to that of 
no magnetic field case is not very clear. It may be possible that the motion of plasma 
in a concave magnetic field (position R) causes Rayleigh-Taylor instability [206] to occur 
resulting in sharp fall in emission intensity. 

5.4.4 Langmuir probe spectra 

Langmuir probe was kept at a distance of 4 mm and at an angle of 45° from the target 
surface where laser beam hit the surface. The output signal was recorded by varying the 
probe voltage from -160 V to +160 V in presence and absence of magnetic field. A triple 
peak structure is observed for the positive probe voltage as shown in Fig.5.27. It was 
suggested [157] that peak 1 is due to photoelectric emission current from the probe caused 
by the UV light emitted by the dense plasma, peak 2 is due to the fast component of 
the plasma and peak 3 is due to the main component of the plasma produced by laser- 
target interaction. It is found that increasing the voltage from 10 Volt onwards intensity 
of peak 1 and 2 remains almost constant whereas intensity of peak 3 increases. The I-V 
characteristics for the case of plasma motion in presence and absence of magnetic fieid 
has been calculated from the Langmuir probe signal using the third peak as shown in 
Fig.5.28. It is observed that in case of magnetic field the ionic part ( the negative voltage 
side) does not saturate unlik e the case of no magnetic field where it saturates. The flow 
of ions are mainl y governed by curvature of the magnetic field which may be responsible 
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Figure 5.27: A typical Langmuir probe signal at a distance of 4 mm from the target and 
at a probe potential of + 10 Volt 



Figure 5.28: I-V characteristics of Langmuir probe signal for the case of plasma in presence 
and absence of magnetic field. 
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Figure 5.29: Square of the probe current is plotted against probe voltage and along with 
a straight line fit. 

for the unusual behavior of ion current not saturating in magnetic field. The square of 
the probe current is plotted with negative probe potential as shown in Fig. 5. 29. for the 
case of magnetic field. The slope of the curve gives the ion density according to the Eqn. 
(2.10). The ion density in case of magnetic field is found to be higher than that of no 
magnetic field case. The estimated density at an angle of 45° to the target normal is in 
agreement with earlier work [100]. The increase in ion density in magnetic field has been 
observed earlier by us using ICCD photography and emission spectroscopy and is due to 
compression of plasma due to magnetic field as discussed before. Eqn. (2.9) is used to 
-calculate electron temperature where ln( \ - L_o ) is plotted against probe voltage for both 
magnetic field and without magnetic field. The electron temperature in case of magnetic 
field at a distance of 4 mm and at an angle of 45° is estimated to be 4.57 eV compared 
to 4.21 eV without magnetic field. Increase in temperature in the presence of magnetic 
field has been observed earlier using emission spectroscopy which is due to ohmic heating 
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caused due to surface current in magnetic field. 


Chapter 6 

Instabilities in a Laser Ablated 
Plasma Expanding in a Magnetic 
Field 


6 . 1 Introduction 

The generation of a dense high-temperature plasma using lasers is of prime impor- 
tance in the field of plasma physics due to its technological applications and rich physics 
un derlying the phenomena. A hot plasma which expands across a magnetic field offers the 
possibilities of investigation of instabilities. The physics of plasma flow across magnetic 
field lines regarding its collimation and stability properties are of particular relevance to 
the solar wind evolution [207], astrophysical jets [208], bipolar flows associated with young 
stellar objects [209], beam plasma heating [10] and cross field injection fueling of tokomak 
plasmas [210]. 

Plasma propagation across magnetic field falls between two extreme limits. In case 
°The work presented in this Chapter is based on : 

1. "Instabilities in laser produced carbon plasma expanding in a nonuniform magnetic field at low pres- 
sure ambient atmosphere" A. Neogi and R. K. Thareja, communicated to Phys. Plasmas. 
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of low /? (particle pressure / magnetic pressure) plasma i.e. in the nondiamagnetic limit, 
magnetic field simply diffuses into whole of the plasma and the plasma moves across B 
field by E x B drift [211] whereas in case of high 3 plasma i.e. in the diamagnetic limit, 
the plasma pressure being higher than magnetic pressure, the plasma "punches” its way 
across the magnetic field giving rise to Rayleigh Taylor (RT) instability as heavy fluid 
(plasma) is decelerated by light fluid (magnetic field) [134, 212]. The Kelvin-Helmhoiz 
(KH) instability is another kind of magnetohydrodynamic instability which arises in a 
stratified fluid if the different layers are in relative motion. Angelo [213] has examined 
the conditions of this instability in a fully ionized plasma immersed in a magnetic field 
and observed KH instability along the magnetic field. Peyser et. a l. [131] designed an 
experiment where plasma is produced by interior illumination of a hollow glass cylinder 
which provides some degree of collimation to the emerging plasma which travels across 
magnetic field giving rise to plasma jet (above 3 k G). Surfaces of these jets undergo edge 
instability, time period of which matches with a modified form of KH instability [139-141.. 

There are three standard methods normally used to determine the stability or the 
properties of the instabilities associated with a given plasma equilibrium as discussed in 
detail in Chapter 1. They are intuitive approach in which the equilibrium is subjected to 
a perturbation that alters the forces acting on the plasma. The second method is based 
on energy considerations where equilibrium state of a plasma will be stable if potential 
energy is a minimum and may be unstable otherwise. A third way to determine the 
stability of a plasma equilibrium is normal mode analysis. This method assumes that if 
the eq uili brium under investigation is perturbed, the linearized plasma equations for the 
time development of the perturbation can be solved, subject to the appropriate boundary 
conditions, assuming a time dependence of the perturbation of the form exp(-i t). 

In this Chapter, we report on propagation of laser produced carbon plasma in an am- 
bient atmosphere in a curved magnetic field (maximum of 3.5 kG). A characteristic delayed 
peak is observed in the temporal profile of CII transition at 426.7 nm and CHE transition 
at 465.0 nm. This delayed peak is not observed in case of only magnetic field or only 
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ambient atmosphere but is observed when both magnetic field and ambient atmosphere 
are present. Moreover, it is observed only in that region of the magnetic field where field 
lines are concave to the moving plasma which is one of the condition for RT instability to 
occur in magnetic field [202]. 

6.2 Rayleigh Taylor Instability 

Rayleigh Taylor instability occurs at the interface of the two fluids of different densities 
when accelerated in a direction perpendicular to their interface; this interface is unstable 
or stable according to whether the acceleration is directed from heavier to the lighter 
medium or vice-versa [202, 206]. In other words Rayleigh Taylor instability occurs when 
the density grows in the direction against the accelerating force. A heavy fluid supported 
against gravity by a light fluid is a common example of Rayleigh Taylor instability. Studies 
of the behavior of small perturbations can be described by [25] 

f = f 0 exp[i(kir - u x t)] ( 6 . 1 ) 

For the motion of plasma (density = p\) in ambient gas (density = po), we get 

u\~ — — koj — — ^ = —ka c a ( 6 . 2 ) 

(Pi + Po) 

Pi - Po 

<*1 = — : — 

Pi + Po 

where a c is the acceleration of the contact boundary. The plasma boundary is stable if us 2 
> 0 whereas if u 2 < 0 the perturbation grows expontially. The growth time of perturbation 
is given by 

r p = (kaai)~% (6.3) 

In plasma, a Rayleigh Taylor instability can occur when a dense plasma (heavy fluid) 
is supported against gravity (light fluid) by the pressure of a magnetic field. Now if the 
plasma is disturbed slightly its motion can deform the magnetic field lines in such a way 
as to produce magnetic forces that tend to amplify the original disturbance. The plasma 



6. Instabilities in a Laser Ablated Plasma 


126 



v ions 

DG 



Mg x Bo 



O B:Z0 





y 



Figure 6.1: (a) Plasma (heavy fluid) supported by magnetic fluid (light fluid) against 
gravity; (b) Ripples at the surface got enlarged due to E x B drift giving rise to Rayleigh 
Taylor instability. 

moving in a curved magnetic field experiences a centrifugal force due to particle motion 
along the curved lines of force which act as an equivalent of "gravitational force" in a 
conventional system discussed above. 

Consider a plasma boundary lying in a Y-Z plane as shown in Fig.6.1(a). Let there is 
a density gradient Ano in the -X direction and a gravitational field g in the X direction. It 
is assumed that 3 is low such that B 0 is uniform throughout. The plasma is supported due 
to magnetic pressure. The ions and electrons initially drift parallel to the Y axis as a result 
of the gravitational force and the magnetic field which is commonly called "gravitational 
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drift" in plasma physics. 


VDG = 


Mg x B 0 

zBl 


(6.4) 


Since the drift velocity is inversely proportional to the charge to mass ratio, the ions and 
electrons drift in opposite directions and the ions have a velocity that is ^ greater than 
electrons. 

If a ripple develops at the interface as a result of random thermal fluctuations, the 
drift v 0 will cause the ripple to grow as shown in Fig.6.1(b). The drift of ions causes a 
change to build up on the sides of the ripple and an electric field develops which changes 
sign on going from crest to trough in the perturbation. This electric field interacts with 


magnetic field giving rise to well known E x B drift, 

Ei x B 0 


vde 


Bl 


(6.5) 


This Ei x Bo drift is upward where the surface has moved upward and downward the 
ripple grows as a result of these properly phased Ei x B 0 drifts. The growth rate of rhe 
instability is calculated using the usual linearized wave analysis. 

The condition for instability to occur is 

(6.6) 


Tin l j2 2 

g— > -k‘v o 
no 4 


where n^ = ^ is the density gradient and k the wave vector. The growth rate is given by 


r n or 
7 S = [~ 9 —] 

rig 


r9^ 

i, = [-] 


£ 


(6.7) 

( 6 . 8 ) 


where s is the density gradient scale length given by ^ . From the condition of instability, 
Eqn.(6.6) the instability occurs when g and ^ have opposite sign. This is similar to the 
statement that light fluid is supporting the heavy fluid. In practical situation the plasma 
is hardly affected by gravity. The effect of magnetic field curvature is similar to g where 
plasma particles flowing along the curved magnetic lines experiences a centrifugal force. 
Thus, if the lines of force are concave to the moving plasma, the centrifugal force and 
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the density gradient are in opposite direction fulfilling the condition of Rayleigh Taylor 
instability. However, if the lines of force are convex to the moving plasma, the plasma tend 
to be stabilized. 

6.3 Results and Discussion 

The plasma expanding in vacuum, low pressure ambient atmosphere, variable magnetic 
field in vacuum and in presence of air is studied by time resolved emission spectroscopy. The 
time resolved emission spectra of carbon plasma were recorded at different distances away 
and parallel to the target surface by moving the monochromator in a plane perpendicular 
to the target surface. The temporal profiles of Cl transition 2p 3 3 D° - 8f F(5/2) at 399.7 
nm, CII transition 3d 2 D - 4f 2 F° at 426.7 nm, CIII transition 3s 3 S - 3p 3 P° at 465.0 
mn and CIV transition at 3s 2 S - 3p 2 P° at 580.1 nm were recorded at different discances 
from the target surface at various pressures of air. Each of these transitions is one of the 
strongest transitions in their respective ionic states. The emission lines are identified using 
the information available in the literature [184]. 

In case of vacuum a single peak structure is observed for all the species Cl, CII, CIII 
and CIV whereas double peak structure is observed for the case of plasma expanding in 
low pressure amb ient atmosphere (5 mTorr, 100 mTorr air) in Cl, CII and CIII transitions 
and has been discussed in detail in Chapter 4. This double peak structure is attributed to 
stratification of the plasma into slow and fast components. In case of variable magnetic 
field in vacuum, the temporal profile of the species show multiple peaks and oscillations 
are observed mainly due to edge effects. It is observed through ICCD photography that 
the plume splits into two lobes which move towards the two poles of the magnet. Each 
lobe contains a fast, intermediate and slow component which gives rise to multiple peak in 
the spectra. A detailed study of the plasma expansion in variable magnetic field in vacuum 
is described in Chapter 5 [129, 178]. 

When air is introduced at low pressure in magnetic field, the multiple peaks and 
oscillations observed in the temporal profile disappear. The nature of the temporal profile 
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Figure 6.2; The temporal profile of Oil transition at 5 mTorr air pressure with and without 
magnetic field. 

is dominated by the gas rather than that by the magnetic field, i.e. the temporal profile in 
5 mTorr/ 100 mTorr air in magnetic field resembles closely the case of 5 mTorr/100 mTorr 
air without magnetic field. The difference caused due to the presence of magnetic field 
is the increase in intensity of the temporal profile, dominance of the intensity of the first 
peak over the second peak for a larger distance and the increase in onset distance of double 
peak compared to no magnetic field case (i. e. only gas). Figure 6.2 shows the temporal 
profile of CIE at 20 mm at 5 mTorr air pressure with and without magnetic field. The 
delayed second peak is of higher intensity in magnetic field at this distance and onwards. 
This is because of ohmic heating at the front edge of the plume which is caused due to 
surface current which arises when plasma moves in a magnetic field. 

An interesting feature is observed when plasma expands in 100 m Torr air in magnetic 
field. The second peak gets suddenly delayed after a certain critical distance in case of CHI 
and CII but not for Cl. In order to investigate the above process the experiment was done 
keeping the target in six different positions in magnetic field at Tl, T2, T3, T4, T5 and 
T6 as shown in Fig.6.3. It is observed that the temporal profile of CIH transitions show 
a sudden increase in delay of the second peak after a certain critical distance as shown in 
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Figure 6.4: Variation of delay of the second peak as a function of distance from the target 
position Tl, T2, T3, T4, T5 and T6 respectively. 

Fig.6.4 for all positions of target except at Tl. Position T2 shows increase in delay after 
16 nun, position T3 after 12 mm, position T4 after 10 mm, position T5 after 12 mm and 
position T6 after 16 mm respectively. It is-- important to note that this delayed second 
peak occurs only at the latter half of the magnet and not in the first half of the magnet. 
Figure 6.3 shows the position of the target and the corresponding point of increased delay, 
e.g. for position T2 it occurs at T2' and similarly for others. Figure 6.5 shows that as 
the target is moved towards increasing side of magnetic field from T2 to T4 where the 
plasma first expands in the convex part of magnetic field and then concave region, the 
onset distance for the delayed second peak decreases. For the position T5 and T6, the 
plasma expands completely in the region where magnetic field lines are concave. As target 
is moved from position T5 to T6 the onset distance for the delayed second peak increases. 
Another interesting feature observed is the appearance of a third peak as one goes from 
position T2 to other positions inside the magnet in the direction of increasing magnetic 
field. This occurs just before the occurence of the delayed second peak. Figure 6.6 shows 
a distinct third peak (P3 in the figure) at 12 mm (dashed line) which disappears at 14 mm 
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Figure 6 . 5 : The onset distance for the occurence of delayed second peak for different 
position of the target. 



Figure 6.6: The dark line represents the temporal profile of CIII transition at 12 mm and 
the dotted line is at 14 mm. PI, P2 and P3 are the first, second and third peak respectively. 
The third peak does not appear at 14 mm. 
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(dotted line) giving rise to a delayed second peak when the target is at position T5; the 
peak eventually disappears as in position T6 where magnetic field is much less than that 
at T5. 

Now we consider various processes responsible for the formation of delayed second 
peak and third peak. They may be due to recombination - three body recombination 
and radiative recombination, charge exchange or magnetohydrodynamic instability. If 
recombination is the cause, the delayed second peak should occur more prominently for 
Cl and CII transition and not for CIII transition. However, it is observed for CIII and 
not for Cl, the recombination cannot be responsible for the formation of delayed second 
peak. In charge exchange process an energetic ion or neutral capturs either an ion or an 
electron and gets transfered to an energetic neutral. Since the delayed second peak; is not 
observed for the neutrals, the probability of charge exchange process is rather low. Let us 
see the possibility of magnetohydrodynamic instabilities for the cause of delayed second 
peak. The growth rate calculations of Rayleigh-Taylor (RT) instability in magnetic field 
show that instability can occur when -g (§^) / no > 0 [202] i.e. when g and density 
gradient scale length of the plasma are in a direction opposite to each other. Since g can 
be used to model the effects of magnetic field curvature, so Rayleigh-Taylor instability can 
occur only when magnetic field fines are concave to the moving plasma. This is exactly 
the situation in our case. The delayed second and third peak occur only at the latter half 
of the magnet where magnetic field fines are concave to the moving plasma. When the 
target is kept to the extreme left (position A) such that plasma mostly travel through the 
first half of the magnet where magnetic field fines are convex to the moving plasma, the 
delayed second peak is not observed. The delayed second peak for CIII ions is formed due 
to two consequitive processes. As RT instability grows, cracks or fissures like structure 
developes at the surface which helps air to diffuse into the plasma giving rise to delayed 
peak as shown schematically in Fig.6.7. Since charged particles play the major role at the 
surface no delayed second peak is observed in neutrals (Cl) . 

An attempt has been made to calculate the growth time of RT instability from the 
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Figure 6.7: A schematic diagram showing the expansion of plasma in magnetic in presence 
of air. The Rayleigh-Taylor instability gives rise to cracks or fissures which allows air ro 
diffuse into the plasma giving rise to delayed second peak. 

intensity plots of ICCD photographs of the plume. The maximum growth rate of instability 
is given by [25]. 

7g = ) l 

Where g e ff is the effective deceleration due to magnetic field and s is the density gradient 
scale length. Figure 6.8 gives ICCD images and the intensity plots at 1000, 1500 and 2000 
ns. It is a plot of total emission intensity of all the ions and neutrals. It is difficult to 
resolve the peaks. The growth time of the instability which is the inverse of growth rare 
given by Eqn.(6.8) can be estimated from the values of g e // and s. The density gradient 
scale length at 1000 ns is estimated to be 5 mm from Fig.6.8 (the length from the plasma 
boundary where density gradient remains constant). Figure 6.9 gives the plot of velocity 
of the plasma front with time calculated from ICCD photographs. g e /f is estimated to be 
0.554 cm fir sec 2 . The growth time of the RT instability which is inverse of growth rate 
comes out to be 0.95 p. sec. 

To conclude, a sudden increase in the delay of the second peak in CIII transition is 
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Figure 6.8: The ICCD images and intensity profile of the plume from ICCD photographs 
at 1000 and 1500 ns and 2000 ns. 
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Figure 6.9: The temporal variation of the velocity of the plume front calculated from ICCD 
photographs. The slope of the curve gives the deceleration in magnetic field. 

observed in case of carbon plasma expanding in nonuniform magnetic field in 100 mTorr 
air by changing the position of the target in magnetic field.lt is observed that the enhanced 
delay occurs in that part of magnetic field where magnetic field lines are concave to the 
moving plasma. Since Rayleigh Taylor instability occurs when plasma expands in concave 
magnetic field, so Rayleigh Taylor instability along with diffusion of air is attributed as 
the cause for the formation of delayed second peak. The growth rate of Rayleigh Taylor 
instability is est ima ted from the intensity and velocity profile of ICCD photographs. 



Chapter 7 
Conclusion 


To conclude, the laser ablated carbon plasma expanding in an ambient gas and 
nonuniform magnetic field has been investigated using emission spectroscopy, fast pho- 
tography and Langmuir probe. Modifications are done in one dimensional Lagrangian 
code MEDUSA to simulate the behavior of the plasma front expanding in ambient at- 
mosphere of He, Ne, Ar and Xe. The simulated result match well with the experimental 
observations. The two dimensional images of the plasma plume expanding in low pressure 
ambient atmosphere exhibit stratification of the plume into slow and fast components. 
The plasma plume expanding in magnetic field shows oscillatory behavior of the plume 
front in agreement with theoretical predictions. Using the time integrated spectra of the 
plasma plume expanding in nonuniform magnetic field it is shown how the gradient and 
curvature of the magnetic field controls the plasma behavior. Multiple peak structure in 
the temporal profiles of various species of the plume is observed in case of plasma ex- 
panding in nonuniform magnetic field. The observed delayed behavior of the temporal 
profile in nonuniform magnetic field in presence of ambient gas is attributed to the onset 
of Rayleigh-Taylor instability in plasma. 

An one dimensional Lagrangian Code MEDUSA which simulates the temporal and 
spatial evolution of laser ablated plasma in vacuum is modified for the plasma expansion 
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in high pressure ambient atmosphere (of He. Ar, Ne and Xe). Rankine-Hugoniot relations 
for the shock wave are used to calculate the boundary conditions of the plasma front. 
Prandtl’s mixing length theory is used to calculate the total shear stress or the viscous drag 
experienced by the plasma boundary. The external gas diffuses into the plasma plume when 
the pressure outside the plasma front is almost equal to the pressure inside the plume which 
gives rise to a diffusive force opposing the motion of the plasma. Finally, the Navier-Stokes 
equation is modified by incorporating the viscous drag term and diffusive force term. The 
simulated results using modified MEDUSA code (with the changed boundary conditions 
and modified Navier-Stokes equation) matches well with the experimental observations. 
The plasma parameters like velocity, temperature, pressure and electron density of the 
plasma front can be simulated for any gas at any higher pressure with the help of the 
modified code. 

The temporal profiles of Cl transition 2p 3 3 D° - 8f F(5/2) at 399.7 nm, CII transition 
3d 2 D - 4f 2 F° at 426.7 nm, CIII transition 3s 3 S - 3p 3 P° at 465.0 nm and CIV transition 
at 3s 2 S - 3p 2 P° at 580.1 nm were recorded at different distances from the target surface 
at various pressures of air (0.1 mTorr, 5 mTorr and 100 mTorr). The intensity of the 
species is observed to increase in presence of ambient gas. The ICCD photographs at 5 
mTorr at various delay times show stratification of plasma plume into two components, 
fast and slow at 600 ns. The velocity of the two components is estimated to be 1.0 x 10 6 
cm/sec and 7.5 x 10 5 cm/sec and hence termed as fast and slow components respectively. 
The ICCD images show that the fast component is bulky at earlier time whereas slow 
component becomes bulky at later time. In vacuum the temporal profiles show a single 
peak structure for all the species. In 5 mTorr and 100 mTorr, double peak structure is 
observed for Cl, CII and CEU transitions whereas no double peak is observed for CTV 
transition. The stratification of the plasma plume is attributed to the diffusion of ionized 
air particles into the plasma plume and collisions with the different charged species of the 
plasma plume. The collided species give rise to slow component whereas uncollided species 
give rise to fast component resulting in stratification. 
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The ICCD photographs of the plasma expanding in a nonuniform magnetic field show 
oscillatory behavior of the plasma front at various delay times (150 and 500 ns) in accor- 
dance with the theoretical predictions. At later time (~ 1000 ns) the plume is observed to 
break into two lobes which move toward the poles of the magnet. The ICCD images were 
used to estimate the plasma parameters. The electron density in presence of magnetic 
field is found to be higher than that of vacuum case and shows oscillations because of 
alternate contraction and expansion of the plasmoid at 150 and 500 ns. The temperature 
of the plasma plume is estimated assuming the plasma plume to behave as an ensemble of 
black body radiators. The temperature in case of magnetic field increases due to Ohmic 
heating caused by the surface current. In presence of ambient gas the temperature rises 
due to enhanced interparticle collisions with the gas particles. The calculated density and 
temperature of the plasma plume match well with the theoretical value. 

The time integrated spectra of the plasma plume propagating in magnetic field were 
recorded at a distance of 6 mm from the target surface keeping the target at three different 
positions in nonuniform magnetic field. It is found that the curvature and the gradient of 
the magnetic field controls greatly the plasma behavior. The plasma expanding in convex 
magnetic field with positive field gradient shows a maximum rise in emission line intensity 
whereas in concave magnetic field with negative field gradient the line intensity decreases 
compared to that of no magnetic field case Furthermore, C 2 molecular Swan band emission, 
d 3 7r s - a 3 7 r u , with A v = -1, 0, 1 is observed to dominate at the former position whereas 
at the latter position it is not observed. The temperature calculated from the ratio of line 
intensity of two lines of same ionization state at a distance of 6 mm from the target is 
higher in former position where the field gradient is positive (11.4 eV) compared to the 
latter where field gradient is negative (4.34 eV). The increase in line intensity is attributed 
to compression of magnetic field lines. C 2 emission in the former position of the target is 
attributed to the movement of electrons due to electromagnetic interaction and that of the 
ions due to curvature gradient drift 

The time resolved spectra of carbon plasma were recorded for all the species Cl, CEt, 
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CIII and CIV in nonuniform magnetic field. It is found that the species at the outermost 
boundary, CII and CIV, show a single peak structure at 4 mm from the target surface 
whereas with increasing distance (~ 10 mm) multiple peaks (six peaks) are observed. 
They are because of two fast, two intermediate and two slow components arising from two 
asymmetric lobes. Multiple peaks and oscillations are observed in the temporal profile of 
CIII and Cl when they are at the outermost boundary (14 mm from the target) where CII 
and CIV disappear. The oscillations are attributed to edge oscillations in magnetic field. 

The temporal profiles of the species show interesting features when the plasma expands 
in combined environment of ambient gas and magnetic field. A sudden delay in second 
peak for CIII and CII species and a triple peak structure is observed but not for Cl. This 
delay is not observed when the plasma expands in either only gas or only magnetic field. 
In order to investigate the delay the target was kept at six different positions in magnetic 
field and it was observed that the delay only occurs when the plasma expands in concave 
magnetic field. This behavior is attributed to the formation of fissures or cracks on account 
of Rayleigh-Taylor instability which occurs when a plasma expands in concave magnetic 
field lines which favors diffusion of air resulting in formation of delayed peak. The growth 
time of the instability is estimated to be 1000 ns from the ICCD photographs. 

Future scope of the work 

It will be an interesting study to simulate the experimental observations of plasma behavior 
in magnetic field using Medusa code. Since J x B and J 2 77 term dominate the behavior of 
plasma motion in magnetic field, so alteration are needed in the energy equation of Medusa 
code. Since the MEDUSA code is one dimensional, two dimensional effects can be studied 
using two dimensional codes like CASTOR. 

The effect of magnetic field curvature on plasma behavior can be studied by varying 
the curvature of magnetic field by changing the distance between the two poles of magnet. 
However, the time resolved spectroscopy of different species at different positions of mag- 
netic field can be carried out. It will augment our knowledge of plasma behavior in curved 
magnetic field. 
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In the present experimental set up the behavior of plasma along the magnetic field 
lines could not be observed because the body of the poles obstruct the line of sight. If the 
plasma motion along the field can be viewed then magnetohydrodynamic instabilities can 
be observed using fast photography. Using those images one can more accurately calculate 
the density gradient scale length to estimate the growth rate of the instability. 

The fast ablated plasma flow into a stationary ambient medium in presence of static 
magnetic field as discussed earlier can be used to study many magnetosheric and astrophys- 
ical processes. At low background gas pressure (P< ITorr) collisionless phenomena such as 
those occuring when solar wind impinges on magnetospheric and ionospheric plasmas and 
with obstacles such as comet, moons, and planets dominate. At high pressure (> 1 Torr) 
collisional interstreaming plasma effects such as blast wave formation are generated which 
behave like supernova shocks. Plasma expanding in magnetic field undergo RT instability. 
An analogous instability was observed in the tail of magnetosphere. The plasma behavior 
in curved magnetic field can be used to study the motion of the plasma in curved magnetic 
field of the earth like generation of ring current in the equatorial plane. Because of cosmic 
abundance of carbon in the interstellar space, spectroscopic study of carbon species is of 
great significance. The knowledge of plasma behavior in different ambient conditions elu- 
cidates our understanding regarding thin film deposition, surgery and many technological 
applications. 
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Appendix A 


c pi. for - used in Chapter 4, 5 and 6 . 



c The following program identifies the peak position, 

c calculates average from different sets of the same 

c distance and finally gives intensity and velocity 

c of different peaks (Peak 1 and Peak2) at different 

c distances. 

c 

c Calculates the peak( 2 only) of given no of points 

c Many sets (2) of input fille 

n 




dimension xl (512 ) , yl (512) 
dimension x2(512), y2(512) 

c COMMENTS 

open (unit =70, access= ' append' , f ile= ' all .out ' ) 

99 format (9X, 'OUTPUT OF 3995 AT 10-3 TORR') 

100 format (10X, 'D ' , 7X, ' II ' , 8 X, ' 12 ' , 10X, 'VI ' , 10X, 'V2 ’ , 

+ 9X, 'PI' , 8 X, ' P2 ’ ) 

write (70,99) 
write (70,100) 
close (unit=70) 

c sc is scale factor. Here taken as 1 

c pi is a parameter for second peak where there is the crest 
c j =0 mean there is only one peak 

c j=l mean there is a second peak, but second peak is smaller 
c j =2 mean second peak is larger than the first 

c DATA 

c 222222222222222222222222222222222222222222222222222222 

d =2 


c 


e=0 . 034 


j=0 

pl=390 


ic=l 


(unit=30 , f ile= ' ol03 . dat ' ) 
(unit=35,file=' 0104.dat') 


CU V L / . . 

11 main (xl,yl/X2 , y2 , d, e, 3 ,pl, sc) 

444444444444444444444444444444444444444444444444 


d=4 

e=0 . 034 
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c 


c 


c 


c 

c 


c 

c 


c . 
5 


j-0 

pl=000 

sc=l 

open (unit=30 , file** ' n905 .dat 1 ) 
open (unit =3 5, f ile= *n906 .dat s ) 
call main (xl , yl , x2 , y2 , d, e, j , pi, sc) 
66666666666666666666666666666666666666666666666666 
d=6 

e=0 .034 
j=0 

pl=340 

sc=l 

open (unit=30,file=' n407.dat' ) 
open (unit =35, file* 'n407.dat' ) 
call main (xl , yl , x2 , y2 , d, e , j , pi , sc) 
88888888888888888888888888888888888888888888888888 
d=8 

e=0 . 017 
j=0 

pl=293 

ac=l 

open (unit=30, file='n409.dat' ) 
open (unit=35, file= 'n410 .dat ' ) 
call main(xl,yl,x2,y2,d,e, j ,pl,sc) 
111111111111111111111111100000000000000000000000000 
d=10 
e=0 . 007 
j=0 

pl=390 

3C= 1 

open (unit=30, file*' n411.dat' ) 
open (unit=35 , file* 'n412 . dat 1 ) 
call main (xl, yl,x2 , y2 , d, e, j ,pl, sc) 
subroutine main(xl,yl,x2,y2,d,e, j ,pl,sc) 


dimension xl (512) ,yl (512) ,x2 (512) ,y2 (512) 
read (30,*) (xl (i) ,yl (i) , i=l, 512) 
read ( 35 , *) (x2 (i) ,y2 (i) , i=l, 512) 

first graph first peak 


call initial (xl,yl,ul,vl) 
ul=ul 


vl=vl 

This above two lines--if there are spikes 

before initial points 

call peak (xl , yl , ul , vl , wl , zl , 0 ) 

call mean (xl , yl , ul , vl , zl , e, wml, zml, 0) 

if (j .NE. 2) goto 5 


tfll=Wl 

Zll=Zl 

tfmll=wml 

Here second peak is termed as second peak 

if ( j .EQ.O ) goto 15 

Lf (j .EQ. 2) goto 10 , 

first graph second peak 1>2 


c 
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200 


c 

c 

c 


75 


80 


write (60,*) d,vel 
write (65, *) d, vein 
close (unit= 50 ) 
close (unit-55) 
close (unit= 60 ) 
close (unit= 65 ) 

format 8 X^ 4 ! ? 7 ' 'i ' 1 

return 

end 


-calculates the initial point 

subroutine initial (x, y, u, v) 
dimension x(512) ,y (512) 
i=0 
i=i + l 

if (i.EQ.250) goto 80 

if ( y (i+1) • EQ .y (i) ) goto 75 

i=i+l 

v=y(i) 

u=x(i) 

goto 80 

continue 

v=v 


u=u 

return 

end 


- calculates the peak (1/2) 

subroutine peak(x,y,u,v,w, z,kp) 
dimension x(512) ,y (512) , s (511) 
real ul 

if (kp.NE.0) goto 90 
do 85 i=l, 512 
Y (i) =y (i) -v 
ul=u 

85 continue 
90 i=0 

95 i=i+l 

s (i) =abs (y (i) ) 
goto no 
105 continue 
i=i+l 

110 if ( kp .NE. 2) goto 115 

if ( x(i) .LT. ul) goto 95 
if ( x(i) .GT. u) goto 140 
goto 120 

115 if (x(i) .LT. u) go to 95 
120 if (y(i) .EQ. y(l)) goto 95 
if (i .GT. 505) go to 125 
if ( abs(s(i)) .GT. abs(y(i+D)) then 
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call peak (xl , yl , pi , vl, wll, zll, j ) 

call mean (xl , yl , pi , vl , zll , e , wmll , zmll , j ) 

goto 15 

c first graph second peak — 2>l 

10 call peak (xl , yl , pi , vl , wl , zl , j ) 

call mean (xl , yl , pi , vl , zl , e , wml , zml , j ) 

15 continue 

write (*,*) wml, wmll 

C2 22222222222222222222222222222222222222222222222222222 

c Second graph starts here 

call initial (x2,y2,u2,v2) 
u2=u2 


v2 =v2 

call peak (x2 , y2 , u2 , v2 , w2 , z2 , 0 ) 

call mean (x2 , y2 , u2 , v2 , z2 , e, wm2 , zm2 , 0) 

if ( j .NE. 2) goto 20 

w22=w2 

z22=z2 

wm22=wm2 

zm22=zm2 

20 if ( j.EQ.O ) goto 30 
if (j -EQ. 2) goto 25 
call peak (x2 , y2 , pi , v2 , w22 , z22 , j ) 
call mean (x2 , y2 , pi , v2 , z22,e,wm22, zm22, j) 
goto 30 

25 call peak (x2 , y2 , pi , v2 , w2 , z2 , j ) 

call mean (x2 , y2 , pi , v2 , z2 , e , wm2 , zm2 , j ) 

30 continue 

Cccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
u= (ul+u2) / 2 
v= (vl+v2 ) / 2 
z= (zl+z2 ) / 2 
w= (wl+w2) /2 
wm= (wml+wm2 ) /2 . 0 
zm= (zml+zm2) /2 . 0 
z=abs (z) *sc 
t=abs (wm-u) 
vel= (d/t) * (1 . 0E+08) 
if ( j.EQ.O ) goto 35 
zll* (zll+z22) /2.0 
wll= (wll+w22) /2.0 
wmll= (wmll+wm22) /2 . 0 
zmll= (zmll+zm22) /2 . 0 
zll=abs (zll) *sc 
tll=abs (wmll-u) 
velll= (d/tll) * (1.0E+08) 
goto 40 

35 zll=0 


40 


velll=0 

close (unit=30) 
close (unit=35) 
open (unit=50, 
open (unit=55, 
open (unit=60, 
open (unit =65, 


access* ' append 1 , 
access- 1 a PP en< ^ ' 
access- ’ append 1 , 
access* 1 append' , 


file='il.out ' ) 
file* 1 i2 .out 1 ) 
file='vl.out ' ) 
file='v2.out ' ) 



Appendix A 


160 



s(i+l)=> s(i) 
else 





a (i+1) = y (i+i) 
endif 





go to 105 




125 

continue 





z=s (i) 





i=0 




130 

i=i+l 





if ( x(i) . LT . u) 

goto 

130 


if (y (i) .NE. z) 

go 

to 

130 

135 

w=x (i) 
goto 150 




140 

continue 
z=s (i) 





i=0 




145 

i=i+l 





if (y(i) . NE . z ) 
w=x (i) 

go 

to 

145 

150 

continue 

return 

end 





c calculates the mean of the peak 

c 

subroutine mean (x, y, u, v, z, e, wm, zm, kp) 
dimension x(512) ,y(512) 
if ( kp .NE. 0) goto 155 
ul=u 

155 continue 
i=0 

ze=abs (z) -e 

jl=0 

sw=0 

sz=0 

160 continue 
i=i+l 

if (i.EQ.500) goto 175 
if (kp .NE. 2) goto 165 
if ( x (i) .LT. ul) goto 160 
if ( x (i) .GT. u) goto 175 
goto 170 

165 if ( x(i) .LT. u) goto 160 

170 if ( y (i) -EQ. y(l) ) goto 160 

if ( abs(y(i) ) .LT. ze) goto 160 
q=x(i) 
sw=sw+q 
r=y (i) 
sz=sz+r 
j l=j 1+1 
goto 160 

175 continue 

wm= (sw) / ( j 1) 
zm= (sz) / ( j 1) 
return 
end 
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c . p2 . for - used in Chapter 5 



c The following program is used in Chapter 5 which 

c prepares data file for temperature calculation 

c using relative intensity of lines of the same 

c ionization state. 

c 

c TWO input fille 

£★★★★★★★*★★★★*★★★★★★***********************★***★** 
dimension a(20),b(20) 
dimension x(20), y(20) 
c open (unit=30 , f ile= ' o319 ' ) 

c * * * * * DATA* *************************************** 
n=5 


c 


c 

c 


c 

c 

c 


5 


6 


7 


8 


20 


n=n-l 

a(i) are the intensity of different lines 

a (1) =109 
a (2 ) =28 . 5 
a (3) =3 
a ( 4 ) =2 
a (5 ) =12 . 50 

b(i) are the corresponding energy value or 

of the upper level in eV. 

b (1) =32.20 
b (2) =34 .28 
b ( 3 ) =50 . 97 
b(4) =40 .57 
b ( 5 ) =4 3 . 25 


read (30,*) (x(i) ,yl (i) / i-l»n) 


do 5 i = 1 , n 

y (i) = ALOG(a(l) /a(i+l) ) 
x (i) = b (i + 1) -b (1) 
continue 
n=n-l 

do 6 i=l,n 

y(4 + i)= ALOG (a (2) /a (i+2) ) 

x (4+i) = b (i+2) -b (2) 

continue 

n=n- 1 

do 7 i=l,n 

y (7 + i) = ALOG(a(3) /a(i+3) ) 
x (7 + i) = b (i+3) -b(3) 
continue 


=n-l 

3 8 i=l,n 

(9+1) = ALOG (a (4) /a (i+4) ) 
(9+i) = b(i+4)-b(4) 

“Sit.50.file-m40.dat') 
= 10 

o 20 i=l,n 

rite ( 50 , * ) x ( - 1 -) • Y 

ontinue 

lose (unit=30) 

lose (unit=50) 


stop 


131098 




